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EXTERNAL AERODYNAMICS OF THE MAGNETOSPHERE* 


By John R. Spreiter, Alberta Y. Alksne, 
and Audrey L. Summers 

Ames Research Center 


SUMMARY 


A comprehensive survey is given of the continuum fluid theory of the 
solar wind and its interaction with the Earth* s magnetic field, and the rela- 
tion between the calculated results and those actually measured in space. A 
unified basis for the entire discussion is provided by the equations of mag- 
netohydrodynamics, augmented by relations from kinetic theory for certain 
small-scale details of the flow. 

While the full complexity of magnetohydrodynamics is required for the 
formulation of the model and the establishment of the proper conditions to 
apply at the magnetosphere boundary, it is shown that the magnetic field 
actually experienced in space is usually sufficiently small that an adequate 
approximation to the solution can be obtained by first solving the simpler 
equations of gasdynamics for the flow and then using the results to calculate 
the deformation of the magnetic field. Numerous specific results are pre- 
sented, including those for conditions in the solar wind, the shape of the 
magnetosphere boundary, and the flow and magnetic fields around the magneto- 
sphere. The calculated results are compared with data obtained by numerous 
spacecraft, including Explorer 12, IMP-1, Pioneers VI and VII, and the Vela 
satellites . 


INTRODUCTION 


The flow of the solar wind past the Earth and its magnetosphere consti- 
tutes a problem similar in many ways to the familiar, but complex problem of 
the external aerodynamics of round-nosed bodies in a supersonic stream. It is 
the purpose of this paper to provide a resume of some of the basic ideas, both 
theoretical and experimental, that support this correspondence, and to present 
a connected account of many of the principal quantitative results that have 
been obtained through application and extension of methods developed origi- 
nally for the aerodynamics of reentry vehicles. In addition, the basic con- 
cepts of the theory of the solar wind itself are outlined in some detail in 
order to illustrate the basic unity of the entire theoretical description of 
the solar wind and its interaction with the geomagnetic field. Finally, 
extensive discussions of the internal consistency of the theory and the degree 
to which the results correspond to conditions observed in space are included 
for the purpose of providing an understanding of the overall reliability of 
the theory. 

^Presented at the Summer Institute for Physics of the Magnetosphere, 
Boston College, Massachusetts, June 19-28, 1967. 



THE SOLAR WIND 


Before embarking on the major topic of this 
solar wind and the Earth *s magnetosphere, let us 
and idealizations of the solar wind itself, many 
figure 1. First of all, it should be recognized 
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Figure 1.- The solar wind - a hypersonic gale in interplanetary space. 


with its nominal velocity of 300 to 800 km/sec, and not the Earth’s 30 km/sec 
motion in its orbit around the Sun, that provides the dominant contribution to 
the relative motion of the Earth and the interplanetary plasma. In many ways, 
we may think of the solar wind as a hypersonic gale that blows through inter- 
planetary space. From an overall point of view, the solar wind may be con- 
sidered to result from forces in the interior of the Sun that are too large to 
be constrained by the feeble pressure that can be imposed by the tenuous cool 
gas of the interstellar medium. As the immediate source of the solar wind, we 
may look to the unsteady conditions in the relatively cool visible surface, 
or photosphere, of the Sun. Here are to be found a wide spectrum of acousti- 
cal, hydromagnet ic, and internal gravity waves that grow tremendously in 
amplitude as they propagate upward into an atmosphere of diminishing density. 
The dissipation of these waves, either directly or indirectly after their cul- 
mination in shock waves, results in a substantial deposition of energy gener- 
ally believed to be sufficient to account for the high temperature of the 
solar corona relative to the photosphere. Since the energy is being 

deposited in the form of heat at a greater rate than can be transported away 

by radiation or thermal conduction, convection in the form of the solar wind 
ensues. The average flow starts gradually with insignificant velocities of 
the order of 1 km/sec or less in the lower corona, and increases steadily with 

increasing distance from the Sun until its ultimate velocity is attained 

beyond a few tens of solar radii. The expanding coronal gas carries with it 
the solar magnetic field, stretching the lines of force outward through the 
solar system and enhancing the strength of the interplanetary field far beyond 
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what it would be in a stationary atmosphere or a vacuum. Representative 
values for the properties of the solar wind in the vicinity of the Earth’s 
orbit are summarized in the box in the center of figure 1. 

At great distances from the Sun, the solar wind must merge somehow into 
the interstellar medium. Although the distribution of interstellar gas, and 
its degree of ionization, is both very nonuniform and uncertainly determined, 
temperatures and densities of the order of those indicated on figure 1 are 
frequently quoted (see, for instance, Kaplan, 1966) . It is generally presumed, 
although on the basis of scanty evidence, that a standing shock wave will 
develop through which the supersonic velocities of the solar wind will decel- 
erate abruptly and become subsonic. Estimates of the location of this termi- 
nal shock wave are very uncertain at this time. It is usually placed in the 
range from 50 to 100 AU. A notable exception, however, is the value of 2 AU 
suggested by Brandt (1962) on the basis of his analysis of the directions of 
comet tails (see Belton, Brandt, and Hodge (1963) and Brandt (1967) for later 
comments) . One aspect of this terminal shock wave that does not appear to 
have been adequately considered in these discussions is that the temperature 
attained by the gas behind the terminal shock wave is far in excess of the 
temperatures quoted above for the interstellar medium. Whether the gas cools 
gradually as it continues to flow subsonically away from the Sun by means of 
some process such as radiation not normally included in the analysis of the 
solar wind, remains at high temperature throughout a spherical shell termi- 
nated at still greater heliocentric distances by a second shock advancing 
through the interstellar medium, or behaves in some other manner altogether 
remains essentially unknown. It may be, on the other hand, that the super- 
sonic solar wind terminates without the intervention of a shock wave in the 
manner suggested recently by Faus (1966) through the action of nonadiabatic 
transfer of kinetic energy to the interstellar medium in a roughly constant 
pressure process that takes place between about 1 to 3.5 AU. 

Although geomagnetic evidence had long before led to the belief that 
solar plasma was projected all the way to the Earth’s orbit and beyond, both 
occasionally from flares and more steadily in localized "M" region beams, the 
modern discussion of the solar wind begins with Biermann’s (1951, 1957) analy- 
sis of the behavior of ionized comet tails. An outstanding photograph of a 
more recent comet, Ikeya-Seki, obtained at the Mojave Tracking Station is 
shown in figure 2. Clearly visible are the head of the comet near the horizon 
at the bottom of the photograph and the extremely elongated tail. Although it 
had become customary to ascribe the antisolar direction of cometary tails to 
radiation pressure of sunlight, Biermann showed this explanation to be inade- 
quate in many respects, and proposed in its place the concept of an approxi- 
mately radial high-speed flow of ionized gas from all parts of the Sun at all 
times. It is not difficult to accept this explanation of a photograph like 
that of figure 2, since the general effect is very remindful of smoke from a 
smudge pot or of a wind sock in a strong and steady wind. 

Although the solar corona is not normally visible because of the over- 
whelming brilliance of the light from the photosphere, it is clearly visible, 
even to the unaided eye, at the time of a total eclipse of the Sun by the 
Moon. As a result, much of what is known about the solar corona has been 
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MOJAVE TRACKING STATION 
GMT 12 h 55-59 m 


Figure 2.- Comet Ikeya-Seki - October 29, 1965. 

learned from observations made during the few minutes duration of such events. 
In more recent years, such observations, particularly of the parts closer to 
the Sun, have been augmented by those taken at other times with the aid of a 
coronagraph, a telescope system fitted with special obscuring equipment 
invented by Lyot in 1937. Figure 3 shows a photograph of the corona taken 
during the eclipse of May 30, 1965, from a NASA Convair 990 jet airliner at 
38,000 feet altitude (Smith, Henderson, and Torrey, 1967). Clearly visible 



Figure 3.- Solar corona during eclipse of May 30, 1965, from NASA CV-990 aircraft at 38,000 feet. 



beyond the disk of the Moon are the polar rays, similar in shape to the mag- 
netic lines of force at the poles of a magnetized sphere. More interesting 
from the present point of view is the large outward extension of the corona 
from the equatorial and middle latitudes. Although these plumes appear to 
extend only a few solar radii in this photograph, other photographs taken dur- 
ing other eclipses have shown rather similar features visible to distances as 
much as five or ten times farther from the Sun. Observations indicate that 
the coronal gas in these plumes is both somewhat hotter and many times denser 
than over the poles. 


In spite of the obvious deviations from spherical symmetry in the solar 
corona, and presumably in the solar wind as well, the theory of the solar wind 
is generally based for reasons of mathematical tractability on the assumption 
of radially symmetric flow. The elements of this theory, which was initially 
presented by Parker (1958) and summarized in more developed form by Parker 
(1963, 1965), Liist (1963, 1965), Scarf (1964), Dessler (1967), and others, is 
outlined in figure 4. Perhaps the most fundamental concept is that of treat- 
ing the solar wind, even as it becomes extremely rarefied far beyond the orbit 
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Figure 4.- Elements of solar wind theory (Parker, 1958). 


of the Earth, as an ordinary continuum gas. The basic equations are thus the 
continuity, momentum, and energy equations of gasdynamics for a radially sym- 
metric outflow of a perfect gas from a centered gravitational field. The 
first two of these are indicated on figure 4, in which p is the density, v 
the velocity, and p the pressure of the solar wind, r is the heliocentric 
distance, Mq = 2x1033 g m i s the mass of the Sun, and G = 6.67xl0" 8 cm 3 /gm sec 2 
is the gravitational constant. Lack of understanding of all contributions, 
such as that of propagating waves, to the energy flux precludes writing a 
reliable energy equation to complete the system. In its place, a variety of 
substitutes have been employed. Among the simplest are the extremes of iso- 
thermal flow within the first few solar radii from the Sun, and adiabatic flow 
beyond. Next in simplicity is the generalization to a polytropic relation 
(p/Po)/(p/Po) a , in which a is an arbitrary constant. The boundary 
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conditions usually employed with these equations are supplied by selecting 
representative values, such as those indicated on figure 1, for the temperature 
and density at the base of the corona, and by requiring that the velocity be 
supersonic, or the temperature or pressure approach zero, at very great dis- 
tances from the Sun. To these must be added the requirement that the transi- 
tion from subsonic to supersonic flow must be achieved with a continuous 
acceleration. If this condition is not observed, a forbidden expansion shock 
wave through which the entropy would decrease would be permitted by the 
governing equations. 

Although the equations for the flow of the solar wind are those of gas- 
dynamics, it has been recognized from the outset that the gas is fully 
ionized, and that the relatively weak interplanetary magnetic field would be 
convected about by the solar wind. Combination of the radially symmetric out- 
flow of the gas with the approximately 27-day rotation period of the feet of 
the field lines anchored in the solar surface leads through this concept, for- 
malized by the relation on the bottom of figure 4 for the direction of the 
magnetic field H in terms of the velocity v of the solar wind and the 
angular velocity of the Sun, to an understanding of the normal spiral 

"hose-angle" configuration of the interplanetary magnetic field. Although the 
configuration of the interplanetary magnetic field is thus primarily deter- 
mined by and secondary to the gas flow, it is of vital importance in guiding 
and deflecting the motion of the relatively rare, but very energetic, par- 
ticles that constitute the galactic and solar cosmic rays. 

Introduction of the polytropic pressure-density law permits integration 
in analytic form of the equations for the solar wind shown on figure 4. 

Parker (1963) has shown that the solution indicates the velocity will go to 
zero near the Sun only if the polytropic index a is less than 3/2. Further- 
more, the solution can only be extended to infinity if the ratio 

R oPo/ GM 0Po = 2R 0 n Q kT 0 / GMqp 0 of the thermal energy to the gravitational poten- 
tial energy at the base of the corona is greater than (a-l)/a. The quantity 
R 0 represents the heliocentric distance to the assumed base of the corona, 
and n represents the number density of the ions. For fully ionized hydrogen, 

n = p/mp in which mp = 1.67><10“ 24 gm is the mass of a proton. Parker also 

showed that R 0 p 0 /GN^)p 0 must not exceed l/2a since no stationary solution 
exists beyond this limit. If the thermal energy should exceed l/2a times the 
gravitational potential energy at the base of the corona, expansion of a more 
explosive character would presumably occur. Within these limits, there still 
remains a substantial family of solutions, each of which involves extended 
heating since a < 5/3. 

Integral curves representing the solution, for a polytropic index a 
within the above bounds, for the variation of the absolute value of the veloc- 
ity v of the solar wind with distance r from the Sun are illustrated in 
figure 5. The singularity in the solution at the point where the two heavy 
integral curves cross occurs where the velocity is equal to the speed of 
sound. It is usually located within a few solar radii of the Sun when repre- 
sentative values are used for the boundary conditions. Various of these 
curves, with either plus or minus sign on v to represent outflow or inflow, 
have been used at one time or another to explain or propose significant 
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FOR I < a < 3/2, 


interplanetary phenomena. 
Earliest is the investigation 
of Bondi (1952) of steady 
spherically symmetric accre- 
tion by a star of a perfect 
interstellar gas considered to 
be at rest at infinity. He 
selected, therefore, the indi- 
cated heavy curve to represent 
an accretion flow that smoothly 
accelerates from subsonic 
speeds at infinity to super- 
sonic speeds nearer the Sun. 
McCrea (1956) later added the 
feature of a standing shock 
wave to terminate the super- 
sonic flow before it reached 
the sensible surface of the 
star. 

In the solar system, how- 
ever, Biermann’s analysis of 
the behavior of comet tails was meanwhile strongly suggesting that the flow 
was not only outward from the Sun, but at high speed. In view of this, and 
also the smallness of the confining pressure that could be maintained by the 
very tenuous interstellar gas, and the even more tenuous, although energetic, 
cosmic-ray gas, Parker considered outflow and selected the indicated heavy 
integral curve to represent a continuous acceleration of the solar gas from 
subsonic speeds near the Sun to supersonic speeds farther away. Although prec- 
edents were numerous in gasdynamics as well as other fields, theoretical objec- 
tions were raised by Chamberlain (1961) to the selection of an integral curve 
that passes through the singular point. He advocated instead selection of one 
of the lower set of curves that go to zero both at the origin and at infinity, 

and never attain sonic velocity at any point. Because the resulting flow is 

about 10 times slower at the Earth T s orbit than that indicated by Parker’s 
solar-wind theory, the results were referred to as Chamberlain’s solar-breeze 
theory. Although it was difficult to determine definitively the relative 
merits of the two proposals at the time, data obtained with many spacecraft 
subsequently have confirmed decisively that the supersonic solution chosen by 
Parker is the appropriate one for at least the portion of the solar system 
within the orbit of Mars, at approximately 1.5 AU. 

Figure 6 shows one of the first results to emerge from Parker’s study of 

the solar wind, the run of velocity with distance from the Sun for a variety 

of assigned coronal temperatures. Although the calculations were made with 
the oversimplified assumption of isothermal flow, the results show two charac- 
teristics that have remained essentially unchanged in all subsequent refine- 
ments of solar-wind theory. They are that the solar wind accomplishes most of 
its acceleration within the first 10 or 20 solar radii, and that velocities of 
the order of several hundred kilometers per second at the orbit of Earth are 
compatible with coronal temperatures of the order of a few million degrees 
Kelvin. 



Figure 5.- Sketch of solutions for variation of 
solar wind velocity with distance from Sun. 
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Figure 6.- Effect of coronal temperature on solar wind velocity - 
isothermal model (Parker, 1958, 1963) . 


Many refinements and elaborations of the theory of the solar wind have 
been made during the years since its first announcement. Although it is not 
the purpose of this paper to present a comprehensive review of the theory of 
the solar wind, a brief outline of some of the salient features of successive 
developments of the theory is presented in figure 7. First of all, we may 
note that nearly all theories of the solar wind are based on the equations of 
a one-component perfect gas. The earliest and simplest version of the theory 
was based on the assumptions of isothermal flow for the first several solar 
radii, and adiabatic flow thereafter to infinity. Since the maintenance of a 
constant temperature in an expanding gas requires a substantial supply of heat 
throughout the isothermal region, this version of the theory may be character- 
ized as having an extended heating region. The region in which heat must be 
supplied to maintain the flow is even more extensive in the versions of the 
theory that employ the polytropic pressure-density relationship because the 
divergence of the energy flux is proportional to 5 - 3a, and is hence positive 
for all a < 5/3. Since the polytropic theory is often employed all the way 
from the Sun to infinity, with either constant or variable a < 5/3, the 
extent of the implied heating region may be very large indeed. 
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Figure 7.~ Successive developments in solar wind theory. 
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Although the damping of waves that are simultaneously tending to grow in 
amplitude as they propagate away from the Sun into a region of diminishing 
density has always remained a possible mechanism for the production of the 
required heat, early estimates suggested the possibility that thermal conduc- 
tion may be the principal process, in addition to convection by the solar wind, 
by which energy is conveyed from the wave-damping region to higher portions of 
the corona. According to this view, it is considered that the shock waves 
that provide the principal heating of the corona are, in effect, absorbed com- 
pletely in a thin heating region in the lower part of the corona, and that the 
thermal conductivity of the hot coronal gases is sufficiently great that the 
heat required at greater distances to drive the solar wind can be supplied by 
heat conduction. Later refinements of this approach have been to add effects 
of viscosity, rotation, and magnetic fields, but these do not usually alter 
the principal features of the flow to any substantial degree. 

It has been customary in the theories of the solar wind that include 
effects of thermal conduction and/or viscosity to use the expressions given by 
Spitzer (1956, 1962) to calculate the transport coefficients. The analysis 
leading to these expressions also indicates that the equiparti tion time 
required for protons and electrons to reach the same temperature tends to be 
of the order of a year when conditions representative of those at the Earth’s 
orbit are used. Using this result, Sturrock and Hartle (1966) have observed 
that one should expect the proton and electron temperatures to be quite differ- 
ent in the solar wind, and proceeded to develop a two-fluid model to allow for 
this effect. Although data from Vela II satellite (Strong, Asbridge, Bame, 
Heckman, and Hundhausen, 1966) show that large differences between the temper- 
atures of the protons and electrons may occasionally occur in the solar wind, 
the temperatures measured most of the time tend to be in better agreement with 
the predictions of the previous single-fluid theories than with those of the 
two-fluid theory. Sturrock and Hartle offer several suggestions to account 
for the apparent heating of the solar-wind protons by a nonthermal energy 
source operative in interplanetary space as far out from the Sun as the 
Earth’s orbit. The possibility they prefer is that protons are heated by sto- 
chastic acceleration following the development of turbulence in the solar wind 
by an unspecified instability mechanism. 


MATHEMATICAL REPRESENTATION OF THE INTERACTION OF THE SOLAR WIND 

AND THE GEOMAGNETIC FIELD 


Just as in the case of the solar wind itself, the fundamental assumption 
upon which this discussion is based is that the average bulk properties of the 
flow can be described adequately by the continuum equations for a perfect gas. 
No attempt will be made here to justify this assumption in detail, except to 
remark that the presence of a weak and irregular magnetic field in the inci- 
dent solar-wind plasma appears to be sufficient to couple the motions of the 
particles even in the absence of collisions, and that numerous data obtained 
in space display outstanding agreement with results calculated on this basis. 
Because of the important role of the magnetic field of the Earth in the inter- 
action with the solar wind, the simplest mathematical model capable of 
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adequate description of many of the dominant features is that provided by the 
standard magnetohydrodynamic equations for steady flow of a dissipationless 
perfect gas (Spreiter, Summers, and Alksne, 1966) . 


The governing differential equations are thus as illustrated in figure 8. 
The symbols p, p, S, and v refer to the density, pressure, entropy, and 
velocity of the gas, H refers to the magnetic field, and y = Cp/c v = (N + 2)/N 
in which Cp and c v are constants representing the specific heats at constant 


DIFFERENTIAL EQUATIONS 

%'P y=o 

A(y-y)y+yp*-(l/47r)£xCURL H 

CURL (Hxv)=0, DIV H = 0 

P/Po 
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{ P/Po ] 



CONSERVATION EQUATIONS FOR DISCONTINUITIES 
[A v n] = 0 

[AV n V + (p + H 2 /87r) fi-H n Ht/47r]=0 
[H n V t -H t V n ]=0,tH n ]=0 
[AV n (h + V 2 /2)+Vn H 2 /47T-H n y.H/47r] = 0 

[Q]-Q,-Q 0 

Figure 8.- Hydromagnetic equations for steady flow of a dissipationless perfect gas. 

pressure and constant volume, and N represents the number of degrees of free- 
dom of the gas particles. For a monatomic gas, N = 3, and y = 5/3. The tem- 
perature T is related to the pressure and density by the equation of state, 
p = pRT/y, in which R = (Cp - c v ) y = 8.315xl0 7 ergs/gm °K is the universal 
gas constant, and y is the mean molecular weight nondimensionalized so that 
y = 16 for atomic oxygen. Thus y = 1/2 for fully ionized hydrogen plasma. 

The electromagnetic parts of these equations are expressed in terms of gaus- 
sian units. Since H = B in such units, these two symbols will be used inter- 
changeably in this paper in order to facilitate further comparisons with the 
original sources. 



Although only first derivatives appear in the system of differential equa- 
tions shown in figure 8, the neglected dissipative terms are described by sec- 
ond derivatives. Their neglect requires that the gradients be small. In 
magnetohydrodynamics, as in gasdynamics, however, compressions tend to 
coalesce and steepen into shock waves of such small thickness that they can be 
considered as virtual discontinuities. In addition, attraction between like 
currents tends to cause current distributions to collapse into thin sheaths, 
across which the magnetic field can be considered in the same sense to be 
nearly discontinuous. Mathematically, continuous solutions of the dissipation- 
less differential equations cease to exist, and the flow is no longer governed 
solely by the differential equations shown on figure 8. Mass, momentum, mag- 
netic flux, and energy must still be conserved, however, and the conservation 
equations provided in figure 8 must hold between quantities on the two sides 
of any such discontinuity . The symbol h in the energy equation represents 
the enthalpy, which is related to the internal energy e = c v T according to 


10 



h =CpT = c y T + p/p. Unit vectors normal and tangential to the discontinuity- 
surface are designated by n and t, and components of v and H in these direc- 
tions are indicated by subscripts n and t. The square brackets are used to 
indicate the difference between the enclosed quantities on the two sides of 
the discontinuity, as in [Q] = Qi - Qo where subscripts 0 and 1 refer to 
conditions on the upstream and downstream sides of the discontinuity. 


Five classes of discontinuities are described by the solutions of the 
conservation equations. As indicated in figure 9, those that lie along stream- 
lines (v n = 0) are called tangential discontinuities or contact discontinuities 
according to whether or not the normal component H n of the magnetic field 
TANGENTIAL 


V n -H n = 0, [v t ]*0, [H t ]*0, [/>]* 0 
[p + H 2 /87r]=0 

CONTACT 

V n =0, H n * 0 , [V] = [H] = [p] = 0, [P]*0 
ROTATIONAL 

V n = ±H n /(47TP) i/2 t [v f ] □ [jj t ] /(47T/>) ,/2 

M s tP] = [v n ] = [v 2 j = [H 2 ] = f Hn ]=0 



SHOCK WAVES, FAST AND SLOW 
V n *0, [P]> 0, [ p ] > 0. [H n ] = 0 

(/ ° v n ) FAST - ^VROT. - ( ^VsLOW 

H, and H 2 (^|*||) THROUGH SH0CK WAVES 


Figure 9.- Hydromagnetic discontinuities. 


vanishes. Discontinuities across which there is flow (v n / 0) are divided 
into three categories called rotational discontinuities, and fast and slow 
shock waves. Some properties that distinguish the various discontinuities are 
summarized on figure 9. 


Of the various possibilities, only the tangential discontinuity has prop- 
erties compatible with those required to describe a boundary surface that 
separates the geomagnetic field and the external flow. As in the classical 
Chapman- Ferraro theory of the magnetosphere boundary based on particle con- 
cepts, the condition H n = 0 holds and requires that there is no connectivity 
between the geomagnetic and interplanetary fields. Although the condition 
v n = 0 prohibits flow across the boundary, arbitrary differences in the den- 
sity and the velocity and magnetic field components parallel to the boundary 
are permitted. It is only required that the sum of the gas pressure p and 
the magnetic pressure H 2 /8tt be the same on the two sides of the 
discontinuity . 


Two important parameters that characterize the flow at any point are the 
Mach number M = v/a and the Alfven Mach number M^ = v/A. The former is the 
ratio of the flow velocity to the speed of sound 

a = (dp/dp) 1 / 2 = (yp/p) 1 /2 = (yRT/m) 1/2 
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The latter is the corresponding ratio to the speed A = (H 2 /4ttp) 1 '^ 2 of a 
rotational or Alfven wave propagating along the direction of the magnetic 
field. Values for a and A for conditions typical of those in the solar 
wind in the vicinity of the Earth's orbit are illustrated in figure 10. 
Although it is possible for weak disturbances to travel as fast as 
(a 2 + A 2 ) 1 / 2 , which they do in a direction normal to the magnetic field, the 

SOUND ALFVEN 

d'(2yRT) ,/2 A=(H 2 /477p) l/2 

R * 8.315 X io 7 erg/gm °K p = m p n 



MAXIMUM PROPAGATION VELOCITY OF A WEAK 
HYDROMAGNETIC WAVE' (Q 2 +A 2 ) 1/2 


Figure 10.- Significant hydromagnetic speeds for fully ionized hydrogen. 


maximum propagation speed is much less than the 250 to 800 km/sec speed of the 
solar wind relative to the Earth, It may be inferred immediately that a 
standing shock wave must exist in the flow around the Earth and its magneto- 
sphere. Of the five classes of discontinuities outlined on figure 9, only the 
fast shock wave can be used to represent the Earth’s bow wave, since the mass 
flow through each of the others is too small. 


DIFFERENTIAL EQUATIONS 


V*? V =0 

rtY-v: )Y + Z p =■ 

(V.V)S=0 , 


P/Po 
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CONSERVATION EQUATIONS 
[A>Vn] = 0 

[AVn V + ( p+ B^ln = 0 
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FLUID MOTION 
I APPROACHES 
f THAT OF 
GASDYNAMICS 


MAGNETIC FIELD 
CURL (IHXV) = 0 
DIV H = 0 


[H n Vt-H t V n ] = 0 
[H n ]=0 


MAGNETIC FIELD MOVES 
WITH THE FLUID 


Figure 11.- Decoupling of equations for large Alfven 
Mach number = (4ttpv 2 /H 2 ) */ 2 . 


An important simplification of 
the magnetohydrodynamic equations 
shown on figure 8 can be made on the 
basis of the large magnitude of the 
Alfven Mach number. This is because 
the order of magnitude of the inertia 
term p (v • V)v in the differential 
equation for the momentum is related 
to that of the magnetic term 
(1/4tt)H x curl H by the square of 
the Alfven Mach number. When the 
latter is large, therefore, the mag- 
netic term can usually be safely 
dropped from the equation, as indi- 
cated in figure 11, with little loss 
of accuracy. Although factors of 2 
appear when similar comparisons are 
made of the magnetic and inertia 
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terms in the conservation equations for the momentum and energy, the general 
conclusion that the magnetic terms may be dropped from these equations when 
the Alfven Mach number is large still holds. An immediate consequence is that 
both the differential equations and the conservation equations defining the 
flow are decoupled from those involving the magnetic field. Thus, all the 
properties of the flow except the magnetic field H can be determined by solv- 
ing the equations of gasdynamics, and the magnetic field can be determined 
subsequently by solving the remaining equations using the values for v 
already determined. Although certain ambiguities reside in the statement, the 
latter equations may be interpreted as indicating that the magnetic field 
lines move with the fluid, the general effect being essentially similar to 
that of a line of smoke or dye released at some instant into a medium flowing 
at speeds high relative to the diffusion rate of the contaminant. In this 
way, one achieves the same decoupling of the calculation of the flow from that 
of the magnetic field as was indicated previously in figure 1 for the solar 
wind itself. The main difference is that the equations for the latter are 
greatly simplified by restriction to radially symmetric flow, although compli- 
cated somewhat by the need to consider the gravitational attraction of the Sun 
and nonadiabatic flow. 


At the magnetosphere boundary, a simplification of great practical util- 
ity can be made by introducing certain approximations to the tangential dis- 
continuity conditions. The basis for these is the fact indicated on figure 12 
that estimates of the magnitudes of the gas pressure p and the magnetic 
pressure H 2 /8tt show that p is much less than H 2 /8tt in the outer magneto- 
sphere, and that p tends to be much larger than H 2 /8tt in the flow around 
the forward part of the magnetosphere. As far as the flow exterior to the 
magnetosphere is concerned, the discontinuities at the magnetosphere boundary 
may thus be approximated satisfactorily by those of the limiting case of a 
tangential discontinuity in which there is a vacuum on one side, and no 
magnetic field on the other side. 


It has also been shown (Spreiter, Summers, and Alksne, 1966) that the gas 
pressure of the flow on the magnetosphere boundary can be approximated 


MAGNETOSPHERE BOUNDARY 
(TANGENTIAL DISCONTINUITY) 
v n = = 0 

[V t ]#0, M*O f [H t ]*0 
<F + h 2/ 87T) mag . =(P + H^7 T) flow 

»P st COS 2 >j/ = K/5 co v2 COS 2 if/ 



div H = 0 
p-- o 


PERMITS CALCULATION 
OF MAGNETOSPHERE 
SHAPE INDEPENDENTLY 
OF FLOW 


MATHEMATICAL FORMULATION OF FREE -BOUNDARY PROBLEM 
FOR SHAPE OF MAGNETOSPHERE 



EQUIVALENT TO OLDER 
CHAPMAN -FERRARO 
PROBLEM BASED ON 
PARTICLE CONCEPTS 


Figure 12.- Decoupling of calculations of magnetosphere shape and surrounding flow. 
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adequately by the simple Newtonian formula p = p<. t cos 2 ip where \p is the 
angle illustrated on figure 12 between the outward normal to the boundary and 
the flow direction of the undisturbed incident solar wind, and p st = Kp^v^ 2 
is the stagnation pressure exerted on the nose of the magnetosphere. The sym- 
bols Poo and v ra in the latter expression refer to the density and velocity of 
the solar wind as it approaches the bow wave. The constant K is usually 
taken to be unity, but yields better accuracy if equated to 0.88 for y = 5/3. 
It is implicit in the use of this expression that cos ip >_ 0, since otherwise 
the surface would be shielded from the flow. The net result of the introduc- 
tion of the Newtonian approximation is that the calculation of the shape of 
the magnetosphere boundary is decoupled from the detailed analysis of the sur- 
rounding flow. Of further significance is the fact that the resulting mathe- 
matical problem defined here on figure 12 is identical to that of the 
classical Chapman- Ferraro theory based on quite different concepts of the 
operative physical processes. 

In this way, the theoretical study of the external aerodynamics of the 
magnetosphere may be separated into three simpler parts without incurring 
undue error. First of all, the shape of the magnetosphere boundary may be 
considered to be the same as indicated by the Chapman -Ferraro theory. Second, 
the flow around the magnetosphere may be determined by solving the equations 
of gasdynamics for flow around the shape given by the solution of the Chapman- 
Ferraro problem. Finally, the magnetic field in the region exterior to the 
magnetosphere may be determined by use of the concept that the field moves 
with the fluid. 


THE SHAPE OF THE MAGNETOSPHERE BOUNDARY 


Chapman and Ferraro formulated their theory for the shape of the magneto- 
sphere boundary and the distortion of the geomagnetic field contained therein 
well over three decades ago (see Chapman, 1963, for a review) , but the result- 
ing free-boundary problem is sufficiently complex that even approximate solu- 
tions were not obtained until well into the satellite era. Although exact 
solutions in analytic form had been obtained earlier for related two- 
dimensional problems suggested by the Chapman -Ferraro problem, much of the 
progress toward the solution of the actual three-dimensional problem is asso- 
ciated with the use of the approximation outlined on figure 13. The approxi- 
mation that the magnetic field at the magnetosphere be proportional to the 
tangential component of the Earth's dipole field was first proposed to be used 
for the entire boundary by Beard (I960) , although Ferraro (1952) had previ- 
ously used the same approximation in a more limited sense to provide the basis 
for a brief discussion of the shape of the equatorial trace of the magneto- 
sphere boundary. Through combination with the other condition that 
HV8 t t = p s t cos 2 ip at the boundary and the expression 

H = - ( a e 3H eq/ r3 ) (6 sin 6 + t2 cos 0) 

for the magnetic dipole field of the Earth, in which a^ = 6.37xl0 8 cm is the 
radius of the Earth and H e q = 0.312 gauss is the intensity of the geomagnetic 
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1. REPLACE H n =0 AT BOUNDARY BY APPROXIMATION THAT H = 2f (H d| ) f 

WHERE (H dj ) IS TANGENTIAL COMPONENT OF THE GEOMAGNETIC 

DIPOLE FIELD AT THE MAGNETOSPHERE BOUNDARY 

2. COMBINE WITH H 2 /8 tt = Kp m V«| COS 2 t TO OBTAIN 

'i±3cos£i / — ! — mV + ' /sin e + . 

R 6 yR Sin 6 d<t>J R 6 y R d<t>) 

(. , . . sin A COS© dR COS <f> dR \ 2 
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. 1/6 

WHERE R = r/D, AND D = Q e (f 2 H£ q /2 7 T V 2 ) 

3. IN EQUATORIAL PLANE, 0=7772, dR/d0 = O 

dR „ fn 6 sin d> cos <b + ,/r6-i \ , , . 

33 - = R I ^ FOR 7T/2 < d> < 37T/2 

d <#> \ R 6 C0S 2 <#>-| / 

4 . IN PRINCIPAL MERIDIAN PLANE, <£ = ± 7 r/ 2 t dR/d</> = 0 

— =(R tan 0)-^^ OR, SOLVING, R COS 0 = 
d0 R 3 ±2 R 3 *! 

Figure 13.- Approximate solution of Chapman-Ferraro problem; dipole axis perpendicular to 

solar wind direction. 

field at the magnetic equator, the partial differential equation shown on fig- 
ure 13 may be derived for the case in which the dipole axis is perpendicular 
to the direction of the incident solar wind. This equation, which was first 
given by Davis and Beard (1962) , is expressed in terms of a spherical coordi- 
nate system r,©,*? with origin at the center of the Earth, polar axis alined 
with the geomagnetic axis, and the angle 0 measured similar to colatitude 
with respect to the geomagnetic pole situated geographically at 78.6° north 
latitude and 70.1° west longitude near Thule, Greenland. The angle cp is 
measured with respect to a line through the origin that is normal to both the 
dipole axis and the free-stream direction, and is equal to tt/ 2 when directed 
in the upstream direction, and 3 tt/ 2 when directed downstream. 

Considerable simplification occurs if attention is confined to either the 
geomagnetic equatorial plane or the principal meridian plane defined by the 
dipole axis and the direction of the incident solar wind. In each case, the 
partial differential equation reduces by reason of symmetry to an ordinary dif 
ferential equation as indicated on figure 13. The equation for the equatorial 
trace of the magnetosphere boundary was integrated numerically by Beard (1960) 
and again somewhat more accurately by Spreiter and Briggs (1961, 1962). The 
equation for the principal meridian trace was derived by Beard (1960) who used 
the solution R = 1 to show the upstream portion of the curve to be circular. 
Spreiter and Briggs (1961, 1962) subsequently considered the entire family of 
integral curves illustrated in figure 14 that are indicated by the complete 
solution written at the bottom of figure 13, in which K is an arbitrary con- 
stant of integration. They showed that only one possibility exists for a con- 
tinuously connected boundary that remains at a finite and nonvanishing 
distance from the origin in the upstream direction and does not violate the 
condition that cos ip >_ 0^ It is composed of part of the circular integral 
curve R = 1 defined by K = 0 in the left plot joined onto part of the 
integral curve on the right plot defined by K = 3/2 2/3 = 1.890. 
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Figure 14.- Integral curves in the principal meridian plane; dipole axis perpendicular to 

solar wind direction. 

The resulting traces of the magnetosphere boundary in these two planes 
are shown in figure 15, together with a plot providing the geocentric distance 
D to the magnetosphere nose defined in figure 13. The latter shows, for 
instance, that D is about 10 Earth radii when the solar wind has a density 
of 10 protons/cm 3 and a velocity of 300 km/sec. 

Coordinates for the remainder of the magnetosphere boundary were obtained 
subsequently by Briggs and Spreiter (1963) , who solved numerically the partial 
differential equation shown in figure 13. The results are illustrated in fig- 
ure 16. The pronounced dent in the boundary in the vicinity of the principal 
meridian plane is associated with the neutral point at which the magnetic 
field just inside the magnetosphere boundary vanishes in the idealized theory. 
Since the boundary condition H 2 /8tt = KpooVco 2 cos 2 ip requires that cos \p 
must vanish at such a point, it follows that the magnetosphere boundary must 
locally be parallel to the free-stream direction of the solar wind. This 



Figure 15.- Traces of magnetosphere boundary in magnetic equatorial plane and meridian plane 
containing dipole axis and solar wind velocity vector; X = 0. 
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Figure 16.- The boundary of the geomagnetic field; X = 0. 

point has further significance because there is only one field line in each 
hemisphere that extends from the Earth to the magnetosphere boundary in the 
idealized theory, and it joins the magnetosphere boundary at the neutral point. 
Further, all field lines in the magnetosphere boundary diverge from one neu- 
tral point and converge to the other. 


The above discussion has been confined to the case in which the dipole 
axis is normal to the direction of the solar wind, or, equivalently, in which 
the angle A between the direction of the solar wind and the geomagnetic 
equatorial plane is zero. Because of the 23.5° tilt of the Earth’s geographic 


or spin axis with respect to the normal 



Figure 17.- Meridian traces of magnetosphere boundary 
for various angles between dipole axis and direc- 
tion of solar wind. 


to the plane of the ecliptic, and the 
11.6° angle between the Earth’s geo- 
graphic and geomagnetic axes, the 
angle A, however, is not always 
small. It is of interest, therefore, 
to determine the resulting effect on 
the shape of the magnetosphere bound- 
ary. The only theoretical results 
that are available on this topic are 
those for the principal meridian 
plane given by Spreiter and Briggs 
(1961, 1962) . 

Figure 17 presents the results 
of their analysis for several rela- 
tive orientations of the dipole axis 
and the solar wind direction. They 
are determined by joining integral 
curves defined by the solution^ 

R cos(0 - AS) + (1/R 2 ) cos 0 = K of 
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the differential equation 

[R 3 cos (6 - AS) ± 2 cos 0]dR/d0 = R 4 sin(0 - AS) + R sin 0 

in a manner similar to that described for A = 0. The symbol K in these 
expressions is again an arbitrary constant of integration and 
S = sin 9 = sin(±7r/2) = ±1. 

The results shown in figure 17 indicate that the magnetosphere boundary, 
as viewed in geomagnetic coordinates, is influenced substantially by the 
direction of the solar wind. It is important to recognize that the effect is 
largely that of rotation into alinement with the direction of the solar wind. 
This point is illustrated in figure 18 by plotting essentially the same 
results as in figure 17 in terms of a coordinate system alined with the direc- 
tion of the solar wind. They are also organized in a different way by 


EQUINOX SOLSTICE 



Figure 18.- Diurnal and seasonal variations of the boundary; coordinates fixed with 

respect to the Sun-Earth line. 


presenting together the results for the extremes of the diurnal variation at 
two significant times of the year, the equinox and the solstice, as indicated. 
It may be seen that the difference between the results is small, and that the 
only feature that is affected to any significant degree is the location of the 
neutral points. These results also indicate that the large open tail always 
extends exactly in the downstream direction, as it must in the exact solution 
of the Chapman- Ferraro problem. Also in accordance with the exact solution is 
the result that the asymptotic width of the magnetosphere tail in the princi- 
pal meridian plane, for any given A, is always exactly twice the distance 
between the two neutral points projected onto the normal to the direction of 
the free stream. 

An evaluation of the accuracy of the results determined through use of 
the approximation outlined in figure 13 can be gained by comparison with the 
higher order approximations to the exact Chapman- Ferraro problem given by 
Midgley and Davis (1963) and Mead and Beard (1964) for the shape of the entire 
magnetosphere boundary for A = 0. Such a comparison is presented in fig- 
ure 19 of the results for the equatorial and principal meridian planes. All 
three sets of results may be seen to be in reasonably good agreement, although 
the higher approximations tend to indicate a slightly larger flaring of the 
equatorial trace of the magnetosphere tail than indicated by the simpler 
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Figure 19.- Comparison with higher order approximations; X = 0. 

calculations of Ferraro (1952), Beard (1960), and Spreiter and Briggs (1961, 
1962) . A similar trend has also been noted in some comparisons with observa- 
tional data. It is difficult to be certain that this is the proper explana- 
tion, however, because, as will be illustrated later herein, the Newtonian 
pressure formula in the Chapman-Ferraro theory tends to underestimate the 
pressure exerted by the solar wind on the flanks and tail of the magnetosphere. 



Figure 20.- Comparison of exact and approximate 
solutions for analogous two-dimensional 
problem. 


Although no higher order approx- 
imations have yet been determined 
with which they may be compared, it 
is not anticipated that the accuracy 
^ of the approximate results for the 
magnetosphere boundary shown in fig- 
ure 17 deteriorates particularly as 
X departs from zero. Support for 
this belief is provided by the com- 
parisons shown on figure 20 between 
the exact solutions given by Zhigulev 
and Romishevskii (1960) , Dungey 
(1961), and Hurley (1961) and the 
approximate solution given by 
Spreiter and Briggs (1961, 1962) for 
the corresponding problem of the 
interaction of the solar wind and a 
two-dimensional magnetic dipole. The 
latter has been determined in a way 


that is entirely analogous to that outlined on figure 13, except that the mag- 
netic dipole field is now given by H = - (a e 2 H 6 q/r 2 ) (6 sin 0 + r cos 0) and 
the length D used to nondimens ionalize the coordinates of the magnetosphere 
boundary is changed as a result to D 2 = a e (f 2 Hiq/2TrKpooVoo 2 ) 1 ^ . Since the 
purpose of this comparison is to illustrate the aegree to which the shape, but 
not particularly the size, of the magnetosphere can be determined by using the 
approximation H = 2f(H c j^)- t , the coefficient f has been set to 0.913 in 
order to achieve a good fit with the exact solution at one point near the nose 
of the magnetosphere. If f had been arbitrarily equated to unity, as it is 
usually done in the three-dimensional applications, the only change is that 
all coordinates of the magnetosphere boundary would be increased by the square 
root of 1/0.913, or by about 4.5 percent. 
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THE FLOW AROUND THE MAGNETOSPHERE 


The preceding discussion has indicated how the exterior magnetosphere 
flow problem can be reduced without undue loss of accuracy to a purely gasdy- 
namic problem of supersonic flow past a given blunt body, the shape of which 
is determined by solving the standard Chapman-Ferraro problem. Although this 
represents an enormous simplification, difficulties associated with the non- 
linear and mixed elliptic-hyperbolic character of the governing partial dif- 
ferential equations are sufficiently great that it remains necessary at the 
present time to approximate the magnetosphere boundary with an axisymmetric 
shape. In the results to be reviewed here, the shape selected is that 
obtained by rotating the equatorial trace of the boundary given by the approx- 
imate solution of Spreiter and Briggs (1961, 1962) illustrated in figures 15 
and 16. Inspection of the results shown in figures 15, 16, 18, and 19 shows 
that this is a reasonable approximation, except possibly where the boundary is 
dented inward around the neutral points. Even this exception may be of less 
importance than might be supposed, because, as will be shown later, these 
dents result from use of the Newtonian pressure approximation, and would dis- 
appear in any internally consistent analysis based on continuum fluid concepts. 


The gasdynamic equations used to calculate the flow around the magneto- 
sphere and the cylindrical coordinate system in which the results are pre- 
sented are summarized in figure 21. The equations are just those indicated on 


DIFFERENTIAL EQUATIONS 


V •/0V = ° 
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CONSERVATION EQUATIONS 
[/>Vn] =0 

[p^l + pn]=o 
[ / ov n (h 2 +v 2 /2)]=o 

Figure 21.- Gasdynamic equations for steady flow of a dissipationless perfect gas. 

figure 11 after dropping the crossed out terms and introducing the assumption 
of axisymmetric flow. Because of the technological importance of this problem 
in the design of reentry vehicles, as well as intrinsic aerodynamic and mathe- 
matical interest, a considerable effort has been expended on this problem for 
more than a decade, and a variety of methods are now available for its numeri- 
cal solution. The calculations of Spreiter, Summers, and Alksne (1966) to be 
reviewed here have been accomplished with the procedures described by Inouye 
and Lomax (1962) founded on the basic method of Van Dyke (1958) and Van Dyke 
and Gordon (1959) as modified by Fuller (1961). As illustrated in figure 21, 



WHERE [ Q ] = Q,-Q 0 
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the method used for the subsonic region near the magnetosphere nose and the 
immediately adjoining portion of the supersonic region is an indirect one in 
which the location of the bow shock wave and the free-stream conditions are 
assumed known and the associated flow field and body shape are found as part 
of the solution. The desired solution for the specified body, the magneto- 
sphere boundary in the present application, is then found by iteration follow- 
ing judicious selection of the initial trial shape for the bow wave based on 
experience with a vast number of cases of aerodynamic interest. A detailed 
discussion of mathematical aspects, such as convergence, stability, and accu- 
racy of the entire method, has been given by Lomax and Inouye (1964) . The 
solution for the remainder of the supersonic region is accomplished as indi- 
cated by the method of characteristics described explicitly for this 
application by Inouye and Lomax (1962) . 



Figure 22.- Streamlines and wave patterns for super 
sonic flow past the magnetosphere; M = 8, 

Y = 5/3. 


Figure 22 shows a plot of the 
results obtained in this way for the 
magnetosphere boundary and the bow 
shock wave, together with several 
streamlines and characteristic or Mach 
lines of the flow for a representative 
set of values for the free-stream Mach 
number Moo and ratio of specific 
heats y. The Mach lines, which cor- 
respond to standing compression or 
expansion waves of infinitesimal ampli 
tude, cross the streamlines at such 
angles that the local velocity compo- 
nent normal to the wave is exactly 
equal to the local sound speed. Mach 
lines thus exist only where the flow 
is supersonic; they are absent from 
the vicinity of the magnetosphere nose 
because the flow there is subsonic. 


Contour maps showing lines of constant density, velocity, and temperature, 
each normalized by dividing by the corresponding quantity in the incident 
stream, are presented in figure 23. As indicated, a single set of contours. 



Figure 23.- Density, velocity, and temperature fields for supersonic flow past the magnetosphere; 

M ro * 8, Y = 5/3. 
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although with different labels, serves for both the velocity and the tempera- 
ture. These results show that the density ratio p/Poo remains near the maxi- 
mum value (y + 1) / (y - 1) = 4 for a strong shock wave in a gas with y = 5/3 
along nearly the entire length of the portion of the bow wave shown. The gas 
undergoes a small additional compression as it approaches the stagnation point 
at the magnetosphere nose and then expands to less than free-stream density as 
it flows around the flanks of the magnetosphere. The velocity remains less 
than in the free stream, however, throughout the same region. Of particular 
interest is the tremendous increase in temperature of the solar wind as it 
passes through the bow wave. If, for example, the temperature of the incident 
solar wind is 10 5 °K, the temperature at the stagnation point is indicated to 
be 2.23xl0 6 °K. Such a value is of the same order of magnitude as the tempera- 
ture possessed by the gas in the solar corona before it is accelerated to the 
high velocities characteristic of the solar wind in the vicinity of the Earth’s 
orbit. This result serves to remind us that the energy content of the solar 
wind is, to a large extent, dependent on the temperature of the corona. 

As noted in figure 11, the deformation of the interplanetary magnetic 
field by the flow around the magnetosphere can be accomplished directly, once 
the flow is determined, either by direct integration of curl(H x v) = 0, 
div H = 0, and [H n ] = 0 or by considering the field lines to drift with the 
fluid. In general, the field lines are spatial curves, but simplicity may be 
achieved at the expense of completeness by restricting attention to the plane 
through the center of the Earth that contains the velocity and magnetic field 
vectors in the incident stream. Results of two such calculations are shown in 
figure 24 for the same conditions as in figures 22 and 23, namely, = 8 and 
Y = 5/3. The field lines are as given previously by Spreiter, Summers, and 
Alksne (1966) , and the contour lines providing the intensity of the field in 
the plot on the left for 90° angle between the directions of the flow and 
field in the incident stream are as given by Alksne (1967) . The contour lines 
in the plot on the right for 45° angle between these two directions have not 
been presented before, but have been calculated in exactly the same way as 



Figure 24.- Magnetic field in plane of free-stream velocity and magnetic field vectors for 
supersonic flow past the magnetosphere; = 8, y = 5/3. 
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those given by Alksne (1967) for 30° and 60°. The corresponding results for 
0° are not shown here, but can be obtained readily from the gasdynamic results, 
since it can be shown that the magnetic field is parallel to the flow and has 
such intensity that the ratio H/H^ is equal to pv/pooVoo at all points. 


PLAN VIEW 
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FIELD LINE 


FIELD LINE 



Figure 25.- Magnetic field for supersonic flow past 
the magnetosphere; = 8, y = 5/3. 


Figure 25 shows the magnetic 
field configuration and intensity 
for the more complex case in which 
the set of field lines considered 
in the incident solar wind is con- 
tained in a plane parallel to the 
plane of symmetry but offset from 
it by a distance D/3. The lines 
of intersection of the offset 
plane with the bow wave and the 
magnetosphere boundary are indi- 
cated by the solid and dashed 
lines, respectively. Correspond- 
ing results for the case in which 
the offset is 2D/3 are shown in 
figure 26. Since the field lines 
are no longer confined to a plane 
behind the bow wave, the results 
for each case are presented in 
the form of a pair of projections. 
It is of interest that the 
extreme crowding together of the 
field lines at the magnetosphere 
nose, so noticeable in figure 24 
for the plane of symmetry, does 
not appear in figures 25 or 26. 
Instead, the field lines appear as 
if they have slipped around the 
nose of the magnetosphere. 

The results shown in fig- 
ures 24 through 26 also indicate 
that the intensity of the mag- 
netic field is usually a few 
times larger near the magneto- 
sphere nose than in interplanetary 
space, and that it may be either 
larger or smaller than in the 
incident stream along the flanks 
of the magnetosphere. They also 

show that the direction of the field in the region between the bow wave and 
the magnetosphere is highly distorted by the flow, and bears little direct 
relation to the direction of the field in the incident solar wind. 




Figure 26.- Magnetic field for supersonic flow past 
the magnetosphere; M ot = 8, y = 5/3. 


The results shown in figures 22 through 26 are all for = 8 and 

y = 5/3. Since both the velocity and the temperature, and hence the speed of 
sound, of the solar wind vary substantially in the course of time, it is of 
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interest to review some of the results given by Spreiter, Summers, and Alksne 
(1966) for other Mach numbers. Figure 27 shows plots of contours of constant 
density ratio p /p for free-stream Mach numbers of 5 and 12. Although dif- 
ferences are clearly evident, they remain sufficiently small to be of only 
secondary importance for most purposes. The corresponding results for the 
velocity and temperature are shown in figure 28. Like the density, the veloc- 
ity field is hardly affected by the change from Mach number 5 to 12. The tem- 
perature fields for the two Mach numbers are substantially different, however, 
with the temperatures for Mach number 12 generally being about five times 
larger than for Mach number 5. Such strong effects on the temperature can be 
understood easily in terms of the expression for T/T^ shown on figure 23. 

M=5 M=I2 

P'Pa> 



0-1 10 
x/D x/D 


Figure 27.- Effect of Mach number on density field for flow past the magnetosphere; y = 5/3. 


M = 5 V/V m M = 12 



Figure 28.- Effect of Mach number on velocity and temperature fields for flow past the 

magnetosphere; y = 5/3. 


24 



Although the results will not be displayed graphically here, the tempera- 
ture contours can also be put into a form relatively independent of by 

normalizing the temperature by dividing by v^ 2 rather than T^. Optionally, 
we may consider the ratio 2c v (T - ToJ/Voo 2 of the change in internal energy 
Cy(T - Too) to the kinetic energy v ao 2 /2 y both per unit mass, or even better, 
perhaps, the ratio 2cp(T - T^/v^ 2 of the change of the enthalpy 
Cp(T - Too) = yc v (T - T^) to the kinetic energy. That this is so may be demon- 
strated by combining the expression (T - TaJ/T*, = [(y - l)/2]]^ x , 2 (l - v 2 /^ 2 ) , 
the relation h = CpT defining the enthalpy per unit mass, and the relation 
Moo = v 00 /(2yRT 00 ) 1 ^ 2 for the free-stream Mach number of fully ionized hydrogen. 
The quantities Cp and c v in these expressions represent the specific heats at 
constant pressure and volume. They are related to the gas constant R accord- 
ing to R = c v (y - l)/2 = Cp(y - l)/2y = 8.315xl0 7 ergs/gm °K. The resulting 
expression is 2cp(T - T co )/v oo 2 = (1 - v 2 /v co 2 ) . Since v 2 /v 00 2 tends to be 
independent of Mach number at large M^, as illustrated in figure 28 for 
Moo = 5 and 12, and T TO in the same region is very much smaller than T, it 
follows that the ratio 2cp(T - T oo )/v oo 2 also tends to be invariant with 
changes in Mach number at large M^. °° 

The conditions in the flow field are also affected by changes in the ratio 
of specific heats y. Earlier calculations, such as those of Spreiter and 
Jones (1963) , have frequently used y = 2, and Dryer and Heckman (1967) have 
recently used y = 1.2. Neither the theoretical reasons nor the observational 
data provide compelling support, in our estimation, for the use of these 
values. We recommend, therefore, that the usual value for a monatomic gas, 
that is, y = 5/3, be used for the present. It is, nevertheless, of value to 
possess an understanding of the effects of substantial variations of both M^ 
and y on the theoretical results. In order to accomplish this in a concise, 
yet still comprehensive, form, figure 29 from Spreiter, Summers, and Alksne 
(1966) has been included to summarize the effects of M^ and y on the stand- 
off distance of the bow wave from the 
magnetosphere nose. These results, col- 
lected from a large number of calcula- 
tions for a wide range of values for 
Mo, and y, show that the ratio of the 
standoff distance A to the geocentric 
distance to the magnetosphere nose D 
is very nearly proportional to the den- 
sity ratio poo/Pb across the nose of 
the bow wave. The latter is, moreover, 
related to M^ and y by the simple 
expression shown at the bottom of fig- 
ure 29. For Mach numbers large com- 
pared with unity, it reduces to 

P m (y- DM® + 2 Poo/Pb = Cy - 1)/Cy + 1) , or simply 1/4 

P b = (y + D m|> for a S as with Y = 5/3. 
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Figure 29.- Bow wave standoff distance for 
supersonic flow past the magnetosphere. 
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MICROSCOPIC INTERPRETATION OF CONTINUUM RESULTS 


In justifying the use of a continuum theory of fluid flow, be it a gasdy- 
namic, magnetohydrodynamic, or more complicated description, to calculate the 
gross average features of the solar wind and its interaction with the geomag- 
netic field, appeal is usually made to possible interaction mechanisms involv- 
ing electric and magnetic fields, and possibly to plasma instabilities, that 
operate over distances that are not only smaller than the ordinary mean free 
paths of the particles, but, more importantly, smaller than the scale of the 
flow feature being examined. When attention is turned from the gross features 
of the flow to the quantities actually measured on board spacecraft with, say, 
a plasma probe, the situation alters significantly because the instruments are 
much smaller than any of the suggested interaction distances. It thus appears 
necessary to return to a microscopic picture based on kinetic theory both to 
interpret certain features of the data and to carry the theoretical calcula- 
tions to the point where the results are directly comparable with observations. 
A complete treatment of the interaction of the solar wind with the geomagnetic 
field based on kinetic theory is, however, an intractable, and probably 
unnecessarily complex, task at the present time. 

One may, instead, follow the procedure employed by Spreiter, Alksne, and 
Abraham-Shrauner (1966) who combined results for the bulk properties of the 

flow provided by fluid models with relations for the velocity distribution 

provided by kinetic theory. The basic additional assumption is that the 
instantaneous microscopic velocity distribution of the ions is Maxwellian in a 
frame moving with the local bulk velocity of the flow. In ordinary gases, 
near standard conditions, the use of a Maxwellian distribution would hardly 
require comment because its applicability is essentially guaranteed by the 
rapid rate of collisions between particles. In the solar wind, the collision 
rate is very low and it is neither obvious nor necessarily correct that the 
velocity distribution is Maxwellian. Nevertheless, the Maxwellian distribution 
is useful at least as a rough guide, and may even be a reasonably good approxi- 
mation to the conditions encountered much of the time in space. The use of a 

Maxwellian velocity distribution is, moreover, consistent with the use of con- 
tinuum equations of magnetohydrodynamics and gasdynamics in the fluid models, 
and with current practice in analyzing and interpreting data from plasma 
probes on some of the recent spacecraft. 

According to the Maxwellian velocity distribution, the number density dn 
of a species of particles of mass m having velocities between w and w + dw 
in a gas of temperature T moving with bulk velocity "w, equivalent to the 
gasdynamic velocity v, is as given in the first expression on figure 30, in 
which k = 1.38xl0~ 16 erg/°K is Boltzmann's constant, n is the total number 
density of the species of particles independent of their velocity, and u, v, 
and w are Cartesian components of w. We proceed to apply this expression to 
the ions, and note that m = 1. 67x10 ~ 24 gm if they are hydrogen ions, or pro- 
tons. The quantity (2kT/m) 1 / 2 has the dimensions of velocity, and is recog- 
nized to be the most probable particle velocity in a reference frame in which 
the bulk velocity w is zero, that is, in a reference frame moving with the 
local bulk velocity. It is proportional to, and less than, the speed of sound 
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Figure 30.- Equations of Maxwellian velocity distribution. 

in accordance with the relations (2kT/m) 1//2 = (p/p) 1 / 2 = a/y 1 ^ 2 = 0.774a 
which may be derived from the expression a = (yp/p) 1 / 2 for the speed of sound 
by replacing the pressure p by the combined pressure 2nkT of the ions and 
electrons, and the density p by the product mn of the mass and number den- 
sity of the ions, the contribution of the electrons to the density being 
negligible . 

It is useful, because of the axial symmetry of the Maxwellian velocity 
distribution about the direction of the bulk velocity vector, to introduce a 
spherical coordinate system in velocity space related to a rectangular 
Cartesian coordinate system in which the w axis is alined with the bulk 
velocity vector w as illustrated in figure 30. The expression for the 
Maxwellian velocity distribution can thus be rewritten in terms of the par- 
ticle speed c = (u 2 + v 2 + w 2 ) 1/2 and the direction of its motion 9 with 
respect to the direction of the bulk velocity vector w in the form shown in 
the second equation on figure 30. Introduction of the dimensionless quantities 
C and W obtained by dividing c and w by (2kT/m) 1/2 leads to the last equa- 
tion shown on figure 30. This expression shows that the fraction dn/n dC dQ, 
of the particles per unit solid angle and unit dimensionless speed depends 
only on the dimensionless speed C and direction 0 with respect to that of 
the local bulk velocity of the particles being counted, and the ratio W of 
the local bulk speed to the most probable speed in a frame moving with the 
local bulk velocity. The quantity W is_moreover directly related to the 
local Mach number M of the flow since W = w/(2kT/m) i ^ 2 = y 1/2 w/a = 1.29M. 

Plots of dn/n dC dft as a function of 0 for several representative 
values of C and W are presented in figure 31. Also included is a plot in 
the upp£r left part of this figure presenting the location of lines of con- 
stant W for the special case in which M^ = 8 and y = 5/3, as in figures 22 
through 26. The velocity distribution having the corresponding value for W 
is thus appropriate at any point along the contour lines. The velocity dis- 
tributions themselves are general, however, and are applicable as well for 
other values for M^ and y, provided that they be applied at the location 
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Figure 31.- Polar plots of Maxwellian proton velocity distributions corresponding to conditions 

in supersonic flow past the magnetosphere. 


where the designated value for W prevails. The plots included in figure 31 
display how the velocity distribution changes from being completely isotropic 
when W is zero, as at the magnetosphere nose, to being highly collimated 
when W = 10, as in the undisturbed solar wind when the free-stream Mach num- 
ber Mo, = (3/5) 1/2 x 10 = 7.75. For fixed C and W, dn/n dC dn is a maximum 
when 0 = 0 and a minimum when 0 = tr, as can be seen by equating 
d(dn/n dC dft)/d0 to zero and investigating the sign of d 2 (dn/n dC dft)/d0 2 
for 0 = it and 0 = 0. 

In figure 32 the velocity distribution dn/n dC dft is illustrated a 
second way, as a function of the dimensionless particle speed C for fixed 
W and 0. The minimum value for dn/n dC dQ for fixed 0 and W is zero, and 
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Figure 32.- Maxwellian proton velocity distributions corresponding to conditions in supersonic 

flow around the magnetosphere. 


occurs when C is either zero or infinite. The maximum value can be found by 
equating d(dn/n dC dQ)/dC to zero and solving for C. It is 

C = [W cos 0 + (W 2 cos 2 0 + 4) 1 ^ 2 ]/2 

For small W cos 0, dn/n dC dft is thus a maximum when C = 1 + (W cos 0)/2. 

For large W cos 0, the maximum value occurs at either C = W cos 0 or 

C = -1/(W cos 0), depending on whether W cos 0 is positive or negative. 

Since there is only one solution for C for given W and 0, dn/n dC dft 

decreases monotonically on both sides of the maximum, and approaches zero as 

C approaches zero and infinity. 

Combining the results of the two preceding paragraphs, we see that the 
absolute maximum value for dn/n dC dft occurs when 0 is zero and 
C = Cm = [W + (W 2 + 4) 1/2 ]/2. Curves illustrating the variation of dn/n dC dft 
with 0 for C = Cj^j are included on figure 31 and indicated by dashed lines. 
Similarly, the smallest_jnaximum value for dn/n dC dft for fixed W occurs 
when 0 = 7T and C = [-W + (W 2 + 4) 1 ^ 2 ]/2. 

The preceding results can be applied to the flow around the magnetosphere 
by inserting into the equations shown on figure 30 the values for density, 
temperature, and velocity calculated using the fluid models. Figure 33 shows 
the results of a specific application for the case illustrated in figures 22 
through 26, namely, = 8 and y = 5/3. The solar wind is considered to con- 
sist entirely of protons and electrons, and to flow in the undisturbed inci- 
dent stream with a velocity Voo = Woo = 4xl0 7 cm/sec. The temperature T^ in 
the free stream is thus 0.90x10 5 °K. Several kinds of results are included on 
this figure. First of all, there are lines representing the magnetosphere 
boundary and the bow shock wave, streamlines, and contour lines of constant 
speed and temperature, all of which correspond to those shown on figures 22 
and 23. The latter are also contour lines of constant local Mach number M 
and W as indicated on the figure. Small diagrams are presented at several 
points of intersection of the streamlines and the contour lines that 



Figure 33.- Maxwellian proton velocity distributions at several points in the flow around the 
magnetosphere; M = 8, y = 5/3, w = 400 km/sec, T = 90,000° K, w = 10.33. 
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illustrate the velocity distribution of the protons. These diagrams are simi- 
lar to those of figure 31 in that they show on a polar plot the variation of 
number density with direction of motion for particles of selected speeds. The 
speeds selected for each of the diagrams include c = c^ = (2kT/m) 1/2 CM so 
that the curve displaying the absolute maximum value for the number density is 
always shown. The values for c for the other curves are simple fractions 
and multiples of c^. The diagrams differ from those of figure 31 in that 
dn/n^ dc dQ is plotted rather than dn/n dC dft in order to illustrate better 
the relation between the velocity distributions at different points. They are 
also oriented on figure 33 so that 0 = 0 in the direction of the local flow 
velocity. These plots illustrate how the velocity distribution is highly col- 
limated in the incident solar wind, completely isotropic at the stagnation 
point at the magnetosphere nose, and changes back to that of a highly colli- 
mated stream along the flanks of the magnetosphere. Except for a limited 
region in the immediate vicinity of the nose, only a small fraction of the 
total number of particles at a given point move in the upstream direction 
(0 = tt) . Throughout the entire flow, however, the speed cy[ corresponding to 
the maximum number density at each of the stations tends to be relatively con- 
stant in spite of substantial changes in the velocity and temperature of the 
gas. 


EVALUATION OF THEORY AND COMPARISON WITH LABORATORY OBSERVATIONS 


It is informative and important to examine the foregoing theoretical 
results for internal consistency and to compare them with observations made in 
the laboratory, as well as in space. One of the first questions that often 
arises in such a discussion is how well does the simple Newtonian pressure 
formula p = p S £ cos 2 provide for the variation of pressure along the mag- 
netosphere boundary. This can be answered partially in the present applica- 



Figure 34.- Comparison of exact and approximate 
pressure distributions on magnetosphere 
boundary . 


tions by comparing the pressure so 
determined with that indicated by the 
detailed gasdynamic calculations. The 
results of such a comparison, origi- 
nally presented by Spreiter, Summers, 
and Alksne (1966) , are shown in fig- 
ure 34 for several values for and y. 
It can be seen that the simple 
Newtonian expression does indeed pro- 
vide a generally good approximation 
over most of the magnetosphere boundary 
of interest in the present studies. 

The agreement is particularly good over 
the portion of the magnetosphere along 
which the flow is subsonic. No signif- 
icant changes are to be expected, there- 
fore, in the shape of the nose portion 
of the bow wave due to the use of a 
more accurate and necessarily more com- 
plicated calculation of the pressure. 

The Newtonian pressure formula does. 
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however, underestimate the pressure somewhat along the flanks of the magneto- 
sphere, indicating that a revised calculation using the gasdynamic pressure 
distribution would lead to a slightly slimmer magnetosphere. The resulting 
change in the magnetosphere shape would be expected to be small, however, 
because of the magnetic pressure against which the gas pressure is balanced 
according to the boundary condition H 2 /8tt = Vst cos 2 is approximately pro- 
portional to the inverse sixth power of the distance from the center of the 
Earth. Whether a more accurate solution of the complete problem would result 
in a magnetosphere that is slightly thinner or slightly more flared remains 
uncertain, however, because of the counterbalancing difference between the 
exact and approximate solution of the Chapman- Ferraro problem illustrated in 
figure 19. A point of equal interest displayed by these results is the lack 
of significant variation of the pressure distribution with changes in Moo and 
y over the range of conditions included in the calculations. 

The pressure at the nose of the magnetosphere provided by the detailed 
calculations agrees, as it must, with the gasdynamic expression 

Pst/Poo v ~ 2 = [Cy + 1)/2] Cy+1)/Cy '°[y - (Y - 1)/2MJ _1/(y_1) /y 

which simplifies for y = 5/3 to p st /p 00 v 00 2 = 0 . 881 [SM^ 2 / (SM^ 2 - 1)] 3 ^ 2 . 

Since the free-stream Mach number M^ is much greater than unity in the pres- 
ent applications, comparison with the Newtonian expression for the pressure 
at the magnetosphere nose, p st = Kp^v^ 2 as indicated in figure 12, reveals 
that K approaches 0.881. The corresponding values for y - 2 and y = 3/2 
are 0.844 and 0.904. Since the distance D to the nose of the magnetosphere 
is proportional to the inverse sixth root of K, it follows immediately that 
the size as well as the shape of the magnetosphere is very insensitive to vari- 
ations of the Mach number and ratio of specific heats of the solar wind gas. 

One of the assumptions fundamental to the attainment of the numerical 
results presented in the preceding pages is that the Alfven Mach number is suf- 
ficiently greater than unity that the equations of magnetohydrodynamics can be 
approximated satisfactorily by those of gasdynamics for the computation of the 
flow field and that the magnetic field can be calculated subsequently by 
application of the concept of the magnetic field lines moving with the fluid 
as outlined on figure 11. Comparison of the values for the Alfven speed indi- 
cated on figure 10 with representative values for the velocity of the solar 
wind assures that this condition is easily met in the incident solar wind. It 
is more difficult to make a similar a priori comparison for the conditions 
behind the bow wave because both the flow velocity and the Alfven velocity 
vary considerably from place to place. Once the velocity, density, and mag- 
netic field intensity have been calculated for any point, however, it is a 
simple matter to evaluate the Alfven Mach number and determine to what extent 
the required criterion is achieved. The results of such an examination are 
presented in figure 35 for two orientations of the interplanetary magnetic 
field for the case of M^ = 8 and y = 5/3. The results are given in the form 
of a ratio of the local Alfven Mach number M^ to the free-stream Alfven Mach 
number M^ in order to take advantage of the simple proportionality between 
the two quantities at any point in the flow field. If M^ = 10, say, the 
results show that M^ is greater than unity everywhere except in a small 
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Figure 35.- Alfven Mach number distribution; = 8, y = 5/3. 

region near the magnetosphere nose. Although some localized departures from 
the calculated results might be anticipated in this region, it would appear 
that the larger scale features of the flow should be quite adequately repre- 
sented by the gasdynamic and frozen- field calculations. The validity of the 
numerical results should be even more assured as increases beyond 10 and 

the critical region in which is less than unity decreases in size. As 

Ma decreases below 10, however, the critical region expands significantly, 
and increasing differences between the approximate numerical results and more 
exact theoretical or observational results should be anticipated. If 
should ever become as small as 2 or 3, the critical region of small local 
Alfven Mach number would become so extensive that substantial departures from 
the present results should be anticipated. Whether M^ can actually attain 
such low values in the solar wind is not known, but the prospect of finding 
an extensive alteration in the basic flow pattern around the Earth should pro- 
vide an incentive to examine observational data for this particular feature. 

That the procedures employed in numerically solving the nonlinear equa- 
tions of gasdynamics for the flow around the magnetosphere are indeed capable 
of yielding results that compare well with actual aerodynamic measurements is 
illustrated in figure 36. Shown in this figure is a shadowgraph photograph 
of a metal model of the magnetosphere in flight at a Mach number of 4.65 
through an argon atmosphere having a ratio of specific heats y of 5/3. The 
dark shadow of the bow wave is clearly visible somewhat upstream of the nose 
of the model. Superposed on the photograph is a dashed white line indicating 
the theoretical location of the bow wave computed for Mb = 4.65 and y = 5/3. 
It may be seen that the agreement between the actual and computed locations of 
the bow wave is virtually perfect. Although the calculations were performed 
for the actual shape of the model, it should be noted that various experimen- 
tal constraints required the shape of the model to deviate in certain particu- 
lars from the shape of the actual magnetosphere. Although these and other 
deviations should have no bearing on the validity of the present comparisons, 
they will be reviewed for the sake of completeness. 
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The principal constraints arise 
from the fact that the experiment 
was conducted by firing the model 
magnetosphere from a .50 caliber 
light-gas gun into the test chamber 
of the Supersonic Free Flight Wind 
Tunnel at Ames Research Center. In 
ordinary use of this facility, 
models are fired upstream through an 
elongated test section of an other- 
wise normal supersonic wind tunnel. 
The working fluid usually employed 
is air, but other gases can be used 
for a more limited range of test 
conditions that can be reached by 
shooting the models into stationary 
gas. Since the relative velocity 
between any given projectile and the 
gas is thus limited by the allowable 
muzzle velocity, the maximum Mach 
number that can be attained depends 
primarily on the speed of sound in 
the gas in the test section. By 
selecting argon as the gas, it is 
possible to obtain a value of 5/3 
for y and a low enough speed of 
sound that Mach numbers as high as 4.5 or 5 can be achieved by firing the pro- 
jectile into stationary gas. Although these conditions are not identical with 
those experienced by the actual magnetosphere past which flows ionized hydro- 
gen at Mach numbers between, say, 5 and 12 or more, they are sufficiently 
close with respect to the two most important parameters of the analysis, 
namely, y and M^, that the results should be useful and instructive. 



Figure 36.- Comparison of calculated and observed 
positions of bow wave of model magnetosphere 
in flight at Mach number 4.65 through argon; 

Y = 5/3. 


Although the actual magnetosphere has an elongated shape because of its 
extended tail, and is not axisymmetric, the model had to be of modest length 
and axisymmetric to permit firing in the gun. Because of these requirements 
and in order to facilitate comparison with theoretical results, the model was 
contoured to the shape described by the equatorial trace of the magnetosphere 
boundary calculated by Spreiter and Hyett (1963) and Spreiter and Summers 
(1963) using a modified version of the Chapman- Ferraro theory in which the 
right-hand side of the boundary condition H 2 /8 tt = p st cos 2 ^ is replaced by 
Ps + Pst cos2 ^ in which p s is a constant. Results were given for many 
different values for the ratio Ps/Pst* that selected for the experiment is 
p s /p s t “ 0.1. The sunward portion of the resulting magnetosphere is very sim- 
ilar, although not identical, to that for p s = 0 used elsewhere in this 
paper. It follows from the basic rules of subsonic and supersonic flow that 
the location of the bow wave and details of the flow around the forward part 
of the magnetosphere are very nearly the same whether p s /p st = 0 or 0.1. The 
value for this parameter has a profound influence on the shape of the rear of 
the magnetosphere, however, and details of the flow may be expected to differ 
substantially in this region. The latter are not visible in figure 36, how- 
ever, because this particular photograph was taken with a focussing 
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shadowgraph system which has the property of portraying the shock wave with 
good resolution, but at the sacrifice of most other details of the flow. Most 
of the scattered isolated details that are, in fact, visible are without aero- 
dynamic or geophysical meaning and are the result of damage to the viewing 
windows caused by impacting models and their fragments in the course of years 
of use of the facility. 

While comparisons such as those shown in figures 34 and 36 and elsewhere 
tend to confirm the applicability and general consistency of the approximations 
introduced to permit the essential solution of the equations representing the 
fluid approach to the interaction of the solar wind and the geomagnetic field, 
an examination of conditions in the vicinity of the neutral points on the mag- 
netosphere boundary described by Spreiter and Summers (1967) shows that cer- 
tain alterations or refinements of the procedures are required. In a sense, 
the changes in conditions to be anticipated near the neutral points are all 
associated with the dual character of the expression p = Kp^v^ 2 cos 2 ip for 
the pressure on the magnetosphere boundary, and the inadequacy, from the fluid 
point of view, of the results obtained therefrom when the angle ip between 
the direction of the undisturbed flow velocity and the outward normal to the 
magnetosphere boundary approaches 90°. 

On the basis of numerous discussions and analyses of Chapman and Ferraro 
and others, the magnetic field lines comprising the magnetosphere boundary 
have long been considered to turn abruptly through 90° at the neutral points 
and extend to the Earth as a single line. It follows from the equations 
div H = 0 and curl H = 0 listed for the interior of the magnetosphere on fig- 
ure 12 that the magnetic field H must vanish at the neutral points. As is 
frequently observed (see, for instance, Dungey, 1958) , this result when com- 
bined with the boundary condition H 2 /8tt = p st cos 2 ip implies that cos ip = 0, 
or \p = 90°. Since it is also required that cos ip > 0 at all other points, 
this result leads to indentations in the magnetosphere boundary in the vicin- 
ity of each neutral point, as can be seen in the meridian plane traces of the 
boundary presented in figures 15 through 20. 

Completely apart from matters concerning the degree to which the condi- 
tion cos ip = 0 is satisfied, the shape of the magnetosphere boundary must 
actually differ significantly from that illustrated in the vicinity of each 
neutral point. Kendall (1962) discussed the situation within the framework of 
free particle flow, and observed that the shape of the boundary immediately 
downstream from the neutral points must be modified because of the additional 
pressure resulting from the second impact of particles that have ricocheted 
off more upstream elements of the surface. The situation is more complex when 
examined from the fluid point of view, but we may commence the discussion by 
noting that Walters (1966) has pointed out that the presence of a concave 
curve or corner on the magnetosphere boundary in a region of supersonic flow 
would result in the formation of a shock wave attached to the magnetosphere 
boundary. 

The shadowgraph in figure 37 illustrates such an attached shock wave in a 
flow with a free-stream Mach number of 5.03 and a ratio of specific heats of 
Y = 5/3 about a simplified magnetosphere shape with a concave corner. This 
photograph was taken in an experiment conducted in the same way as described 
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for figure 36. The results differ 
because the model, while still axisym- 
metric, was contoured to correspond to 
the principal meridian, rather than the 
equatorial, trace of the magnetosphere 
boundary . 

Although an attached shock wave is 
thus to be expected in supersonic flow 
around a solid obstacle with indenta- 
tions like those of the theoretical 
magnetosphere shapes illustrated in 
figures 15 through 20, it is not what 
is to be expected in high Mach number 
flow around the actual magnetosphere. 
The reason is that the magnetosphere 
is not a solid obstacle of fixed geometry, but a surface that must adjust its 
shape in such a way that the balance between gas and magnetic pressure be 
maintained in accordance with the approximate relation (H 2 /87r) ma g = Pfi ow 
given on figure 12 for a tangential discontinuity. Now the presence of an 
attached shock wave requires that there be a discontinuous change in the gas 
pressure on the magnetosphere boundary at the foot of the shock wave. Since 
div H = 0, curl H = 0, and H n = 0 taken together imply that H is continu- 
ous, there can, however, be no balancing discontinuity in the magnetic 
pressure H 2 /8 tt. 

Additional difficulties occur if the attached shock wave is oblique to 
the magnetosphere boundary, as it appears to be in figure 37, because the 
flow, and therefore the magnetosphere boundary as well, would be required to 
turn a sharp corner. A curl-free and divergence-free magnetic field would be 
either zero or infinite at the apex, depending on whether the corner were con- 
cave or convex toward the magnetosphere interior. Correspondingly, the mag- 
netic pressure would be either zero or infinite at the corner, and a second 
condition would be imposed that would be impossible to achieve with a gas 
pressure. A further difficulty with the fluid extension of the conventional 
theory of the magnetosphere boundary is that the latter indicates that 
H = 0 at the neutral points. The boundary condition expressing the pressure 
balance would therefore require that the pressure p vanish at these points, 
but this is also impossible in a gas. A step toward the resolution of this 
contradiction is suggested by considering a magnetosphere boundary model in 
which the Newtonian pressure formula p = Kp^Voo 2 cos 2 ip is replaced by 
p = p s + KpooVoo 2 cos 2 ip where cos >_ 0 and by p = p s where cos ip <_ 0 . 

As described previously (Spreiter and Hyett, 1963; Spreiter and Summers, 1963; 
Slutz and Winkelman, 1964) , the magnetosphere boundary defined by this model 
contains a sharp cusp directed toward the Earth at each neutral point for all 
values for ps/Kp^Voo 2 = Ps/Pst except zero. Although this result is consis- 
tent with the requirements of the free-particle approach, and has been pro- 
posed again more recently by Grad and Hu (1966) to be the appropriate 
configuration to consider for this region, it is not an admissible solution 
from the fluid point of view because the flow cannot follow the cusped 
contour without developing a shock wave. 



Figure 37.- Focussed shadowgraph of model of 
magnetosphere in flight at Mach number 5.03 
through argon; y = 5/3. 
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Spreiter and Summers (1967) have observed that the difficulties may be 
resolved if the boundary of the flowing solar plasma is considered to separate 
from the boundary of the geomagnetic field so as to form an enclosed region of 
finite extent filled with essentially stationary gas at constant pressure p c , 
as illustrated in figure 38. Since the tangential velocity components are 
different on the two sides of the free streamline that divides the flowing 




Figure 38.- Cusp-shaped regions of hot plasma in the vicinity of neutral points. 

solar plasma from the stationary trapped plasma, and the normal velocity com- 
ponent is zero, conditions on the two sides of the free streamline are con- 
sidered to be in accordance with those of a tangential hydromagnetic discon- 
tinuity outlined on figure 9. Arbitrary changes in density and tangential 
components of the velocity and magnetic field are allowed across such a sur- 
face, but the normal component of the magnetic field must be zero, and the 
sum of the gas pressure p and the magnetic pressure H 2 /8 tt must be the same 
on both sides of the surface. Although the magnetic field is considered to be 
sufficiently strong to provide the necessary collective interaction between 
the individual particles of the rarefied plasma to enable it to act like a 
fluid, it is nevertheless considered to be sufficiently weak that the magnetic 
pressure is otherwise negligibly small compared with the gas pressure on both 
sides of the discontinuity. With previous theory retained as much as possible, 
it follows that p c = p s ^ cos 2 ipf where is the angle ip associated with 

the free streamline, and the magnetic field is such that H c 2 /8tt = p c = const 
and H n = 0 along the portion of the magnetosphere boundary adjacent to the 
trapped plasma. Over the remainder of the magnetosphere boundary, the condi- 
tions specified on figure 12 are to be applied in the usual manner. 

Although the exact solution of this problem has not been determined, an 
approximate procedure has been described by Spreiter and Summers (1967) that 
leads to the results for the shape of the principal meridian trace of the mag- 
netosphere boundary illustrated in figure 38. As can be seen, the streamline 
bounding the flowing plasma coincides with the magnetosphere boundary from 
the stagnation point at the nose to a separation point, situated somewhat 
upstream of the neutral point. The flow then detaches from the boundary of 
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the geomagnetic field and follows a free streamline of constant pressure until 
it reattaches to the magnetosphere boundary at a point somewhat downstream of 
the neutral point. The relation p c = p $1 _ cos 2 \\> = const indicates that the 
free streamline must be a straight line m the principal meridian plane. It 
must also be tangential to the magnetosphere boundary at both the detachment 
and reattachment points since otherwise unwanted discontinuities or singulari- 
ties in the pressure would occur that could not be balanced satisfactorily by 
the magnetic pressure within the magnetosphere. Although the approximate 
results shown in figure 38 do not display the feature adequately, the portion 
of the magnetosphere boundary adjacent to the stationary plasma must be suf- 
ficiently sharply cusped that | H | is constant and different from zero at all 
points . 

Since the gas trapped in the cusped regions must have come originally 
from the solar wind, its temperature T c can be estimated to be the same as 
the stagnation temperature of the solar wind. The latter is given in terms of 
the temperature T^ and Mach number of the incident solar wind by the 

expression T c = T^fl + (y - l)M OQ 2 /2] obtained from the more general expres- 
sion given on figure 23 by setting v/v^ to zero. The temperature of the 
trapped gas is thus independent of the relative orientation of the dipole axis 
and the free-stream direction of the solar wind, and is equal, for example, to 
22.3^ when = 8 and y = 5/3. The number density n c of the trapped pro- 
tons can be estimated by combining the above expression for T c with the 

relation p c = 2n c kT c = p^v^ 2 cos 2 in the derivation of which the coeffi- 

cient K has been equated to unity for simplicity. The resulting expression 
for the ratio n c /noo of the number density of trapped protons to that of the 
incident solar wind is n c /noo = (2 yM ^ cos 2 ipf) / [2 + (y - 1)M 00 2 ] . For flow 

with = 8 and y = 5/3, Spreiter and Summers (1967) find that = 53° and 

n c/ n oo = 1-^3 when the dipole axis is normal to the direction of the solar 
wind. When the direction of the solar wind is inclined 34. 5°, the correspond- 
ing values are 33° and 3.34 in the summer hemisphere and 73° and 0.39 in the 
winter hemisphere. 

Each trapping region may be recognized to have essentially the same shape 
as the axisymmetric cusped geometry configurations considered in studies of 
plasma containment for production of thermonuclear energy. The major advan- 
tage of the cusped configuration is its stability. One of its chief disadvan- 
tages is the large rate of loss of particles. Since a large loss rate could 
have important geophysical consequences, both as a mechanism for injecting 
particles into the magnetosphere in general, and as a source of particles that 
spiral down along the field lines from the neutral points to localized regions 
in the Earth* s atmosphere, it is of interest to inquire into the loss rate 
that might be expected. This has been done by Spreiter and Summers (1967) 
using a basic expression given by Grad (1963) for the proton loss rate u at 
an axisymmetric cusp for an isotropic distribution of particles of a single 
speed V +c and density n Q . It is that u = 45 (m_/m ) 1 / 2 n c V +c 7iX 2 /4 in 
which m_ = 9.11xl0" 28 gm and m + = 1.67xl0~ 24 gm represent the mass of an 
electron and a proton, and A = m + V +c c/eH c , in which e = 4.8xlO” 10 esu is the 
charge on an electron and c = 3x1010 cm/sec is the speed of light, represents 
the gyro radius of a proton moving with speed V +c in a magnetic field of 
intensity H c . Grossmann (1966) has recently reexamined Grad’s problem in 
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greater detail and arrived at an expression for u that differs somewhat in 
numerical coefficient but is otherwise the same as that presented above. The 
corresponding Maxwell averaged loss rate <u> follows directly upon replace- 
ment of n c V£ c by 3nkT/m+ in the expression for u. Insertion of the appro- 
priate relations for conditions in the trapped plasma leads to the expression 
<u> = 2. 8x10 14 [(y - l)/y] ^ 2y 00 - For representative values of 
Voo = 5x 10 7 cm/sec and y = 5/3, this expression indicates a loss rate of 
8.9xl0 21 protons/sec from each cusp. This rate is, moreover, independent of 
the number density of the solar wind and the relative orientation of the dipole 
axis and the direction of the solar wind. The area A = 45ttA 2 over which the 
particle flux is distributed near the tip of the cusp is found to be given by 
A = 2.8x10 16 (y - l)/(Ynoo cos 2 ^f) for free-stream Mach numbers much greater 
than unity. The flux density of protons per unit area <u>/A out of the 
vicinity of the tip of the cusp is thus proportional to the solar wind flux 
n^Voo and exhibits substantial seasonal variations in accordance with the 
values for cos 2 ij;f. With <u>/A normalized to unity when the dipole axis 
is perpendicular to the direction of the solar wind, this result when combined 
with the values for given above shows that the flux density through the 

summer cusp is approximately 2 and that through the winter cusp is approxi- 
mately 1/4.5 when the direction of the solar wind is inclined 34.5° to the 
plane normal to the dipole axis. Since, in addition, the density of the solar 
wind may vary by a factor of 10 or more and the velocity by a factor of 3 or 
more in time intervals as short as a few minutes, changes in the values of 
<u>/A by a factor of 30 or more are to be anticipated in correspondingly short 
time intervals. The total flux of protons from each cusp would vary much 
less, however, because it depends only on v^. 

The possible existence of particle fluxes into the magnetosphere near the 
neutral points appears to be supported by laboratory experiments designed to 
simulate the interaction of the solar wind and the geomagnetic field. Results 
obtained by independent groups in Canada (Osborne, Bachynski, and Gore, 1964), 
Japan (Kawashima and Fukushima, 1964), and USA (Cladis, Miller, and Baskett, 
1964; Waniek and Kasai, 1966) have all displayed injection of plasma into the 
magnetosphere from cusp-shaped regions extending toward the Earth from the 
vicinity of the theoretical neutral points. 

Ground-based observational evidence that such phenomena may actually 
occur is provided by magnetometer records from the polar regions. A consider- 
able body of data gathered during the IGY (July 1957-December 1958) and the 
Second Polar Year (August 1932-August 1933) has been assembled and analyzed 
by Fukushima (1962) in various ways according to the time of day and year, 
and the degree of disturbance of the worldwide geomagnetic field. The results 
show that the polar regions are constantly disturbed, even at periods when the 
planetary geomagnetic field is very quiet. Two of Fukushima T s many summary 
plots are reproduced in figure 39. They show the latitudinal distribution of 
the average geomagnetic disturbance near the noon-midnight geomagnetic merid- 
ians for the case of Kp = Oo, that is, for the very quietest times as indi- 
cated by the planetary geomagnetic index, at two selected sets of times 
extending over the entire period of the IGY. The results show that the magni- 
tude of the polar-cap disturbance is of the order of a few or several tens of 
gammas, that both polar caps are equally disturbed when the geomagnetic dipole 
axis is p erpendicular to the Sun-Earth line (A = 0) , and that the disturbances 
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Figure 39.- Polar magnetic disturbance on very quiet days (Kp = Oo) during I GY (Fukushima, 1962). 

are a maximum at summer noon and a minimum at winter midnight. These results, 
as well as the location of the maximum disturbances, appear to be consistent 
with what might be anticipated on the basis of direct invasion of charged par- 
ticles or propagation of hydromagnetic waves into the polar regions along the 
fields from the neutral points on the magnetosphere boundary. 

The results described above are also consistent with several inferences 
drawn by Lebeau (1965) from analysis of diurnal magnetic activity at a number 
of stations at very high latitudes in the arctic and antarctic. An important 
conclusion of his studies is that stations within a large portion of the polar 
cap display a diurnal variation of magnetic activity that peaks strongly at a 
time midway between local noon and magnetic noon. This is explained in terms 
of a model in which the intensity of the magnetic activity is determined by 
the combined action of two separate factors, an ’’excitation factor” which 
tends to create a maximum at magnetic noon and a ’’modulation factor” which is 

considered to result from solar 
influences on the conductivity of the 
E-region of the ionosphere. The for- 
mer is presumed to be the result of 
particle precipitation, and displays 
many properties that are consistent, 
perhaps even quantitatively, with the 
theoretical results described above. 
From figure 40 reproduced from Lebeau 
(1965) , it may be seen, for instance, 
that the seasonal variation of the 
amplitude of magnetic activity at 
Godhavn in the arctic is almost per- 
fectly out of phase with that at 
Dumont d’Urville in the antarctic 
with the maximum activity occurring 
at each station during the months of 
local summer. Values between about 


controlled by local time. The latter is 
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Figure 40.- Amplitude of magnetic activity for 
an arctic station, Godhavn, and an antarctic 
station, Dumont d’Urville (Lebeau, 1965) . 
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3 to 5 for the enhancement of the amplitude of the activity observed in mid- 
summer months over that at the same station in the midwinter months are compat- 
ible with the seasonal variations of the proton flux density calculated by 
Spreiter and Summers (1967) . Although the latter have not been evaluated 
explicitly for Godhavn or Dumont d'Urville, the seasonal variation should be 
somewhat smaller than the value of 9 quoted above because geometrical factors 
tend to prevent observation at a single station of the maximum value of the 
excitation factor in the entire summer polar cap and the minimum value in the 
entire winter polar cap. The amplitude of the magnetic variation at Godhavn 
may also be seen to vary recognizably with the solar cycle, but with a range 
that is substantially less than the seasonal variation. As an example, mean 
values for July 1962 were only about 70 percent of those for July 1959. Such 
a change would correspond, assuming a linear relation between the amplitude of 
the magnetic activity and the injection rate, to a value for the mean flux 
density PooVoo of the solar wind about 70 percent as large in 1962 as in 1959. 
Although direct measurements of the density and velocity of the solar wind 
were not made in 1959, this result is reasonable in terms of known variations 
in the solar wind and its interaction with the Earth's magnetic field. 


COMPARISON WITH OBSERVATIONS IN SPACE 


The final test of any theory of the interaction between the solar wind 
and the Earth must be based on comparisons with observations actually made in 
space. Such comparisons have already been made by several different groups 
on several levels of detail. These range from simple comparisons of gross 
features such as the existence and location of the magnetosphere boundary and 
the bow shock wave to more detailed comparisons of features such as the veloc- 
ity, density, and temperature of the flowing plasma and the intensity and 
direction of the magnetic field. It is inevitable that even finer details 
such as perhaps the precise form of the velocity distribution of the particles 
as well as temporal fluctuations and changes of the bulk quantities will be 
compared as further development of the theory, instrumentation, and data 
processing techniques continues. 

Many of the earlier comparisons involving the bow wave made use of the 
theoretical results of Spreiter and Jones (1963) shown here in figure 41. 
Included in this figure are the equatorial trace of the magnetosphere boundary 
for X = 0 given previously in figures 15 and 16 together with the correspond- 
ing shape of the bow wave for axisymmetric flow of a perfect gas approaching 
the Eart’ at a Mach number of 8.71. Only the position of the bow wave for 
y = 2 w. given originally, but the corresponding results for y = 5/3 were 
made available shortly later in IG Bulletin 84 (1964). Since it was an 
important part of the original purpose in making these calculations to deter- 
mine the expected location in space of the bow wave and the magnetosphere 
boundary, the results were given in dimensional, rather than in nondimensional , 
form. Values for the number density n^ of the protons and their bulk veloc- 
ity v TO in the incident solar wind were selected as indicated on the figure, 
although the results would be identical for any other combination of values 
that leads to the same product n Qo v oo 2 . They could also be changed simply to 
correspond to any other value for this product by altering all linear 
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Figure 41.- Calculated shape of shock wave and suggested nature of conditions in the transition 
region; n^ = 2.5 protons/cm 3 , = 600 km/sec, B i = 5 gammas. 

dimensions in proportion to the ratio of the inverse sixth root of n^v*, 2 . 

The value of 8.71 for the Mach number was selected so as to be equal to the 
Alfven Mach number corresponding to the values for the density, velocity, and 
magnetic field indicated on the figure. As is evident from the discussion 
associated with figure 11, there is, however, no need to equate the Mach num- 
ber to the Alfven Mach number. A large value for the latter merely insures 
that the equations of magnetohydrodynamics can be approximated by the simpler, 
though still complicated, equations of gasdynamics. A value of 8.71 for the 
Mach number is, nevertheless, both reasonable and sufficiently great that the 
location of the bow wave does not depend critically on the Mach number. 

Also retained on figure 41 are small insert plots suggested to be indica- 
tive of possible appearances of the magnetic field in the region between the 
bow wave and the magnetosphere. These were reproduced from results calculated 
by Auer, Hurwitz, and Kilb (1961, 1962) for the detailed magnetic structure of 
the transition region between a steady magnetic field and an adjacent plasma 
containing a smaller magnetic field. Results of similar calculations for a 
much wider range of conditions have been given subsequently by Rossow (1965, 
1967) and Jones and Rossow (1965) . Although all of these calculations are for 
one- dimensional unsteady flow and the geophysical application is considered 
here as a three-dimensional steady flow problem, the conditions at any given 
instant after the one-dimensional interaction zone has become of sufficient 
width should be approximately equivalent to those actually encountered in 
space between the magnetosphere boundary and the bow wave. These results 
illustrate that a value of 2 for the Alfven Mach number based on the velocity 
component normal to the shock wave is critical, and that conditions associated 
with values on either side of 2 are strikingly different. If M^ < 2, as 
occurs well back along the shock wave in general, and possibly near the nose 
if the interplanetary magnetic field is sufficiently strong, the magnetic 
fields calculated by Auer, Hurwitz, and Kilb (1961, 1962) are oscillatory and 
orderly. A pair of examples for parallel and antiparallel geomagnetic and 
interplanetary fields is shown in the small insert figures for M^ = 1.54. 
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Although the magnetic field variations in the transition region are similar in 
form for the two cases, the results are distinguishable because of the change 
in direction of the field at the magnetosphere interface. Particularly signif- 
icant is the fact that the absolute value of the magnetic field undergoes no 
significant change at the boundary of the magnetosphere. If > 2, the mag- 
netic field profile at any instant is very irregular, as is shown in the 
remaining pair of small inserts in figure 41. They are for an example in 
which Ma = 5.8 and the interplanetary field is either parallel or antiparallel 
to the geomagnetic field. Of particular significance is the fact that the 
irregular field calculated for the transition region retains the same direc- 
tion as the interplanetary field and has an average value that may be much 
less than the intensity of the field at the boundary of the magnetosphere. 


Figure 42 shows some observational results obtained by Cahill and 
Amazeen (1963) with Explorer 12 that bear close resemblance to the results of 
Auer, Hurwitz, and Kilb (1962) for > 2. The continuous line labeled 
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Figure 42.- Explorer XII data resembling calculated results for antiparallel fields with M > 2 

(Cahill and Amazeen, 1963). A 


|B di | represents the values for the intensity of the geomagnetic field calcu- 
lated by harmonic expansion of the surface values according to the analysis of 
Finch and Leaton (1957). Because of the relatively large values for r/a e , 
the calculated intensities are very nearly those provided by the geomagnetic 
dipole term alone, since the other terms are not only smaller at the Earth f s 
surface but diminish more rapidly with increasing altitude. The data points 
and broken curves represent the measured values for the intensity |b| and 
the direction of the actual magnetic field. The angle a lies between the 
direction of the magnetic field and the spin axis of the spacecraft. The 
angle ip lies between the plane containing the magnetic field and the spin 
axis and the plane containing the Sun direction and the spin axis. Each data 
point shown represents the average of 20 individual measurements made at the 
rate of approximately 3 per second, and only a small fraction of even these 
values is included on the plot. Moreover, each component of the individual 
measurements contains an uncertainty of ±12y resulting from digitization of 
the data. The boundary of the magnetosphere, plainly visible at about 
8.2 Earth radii, was traversed by Explorer 12 at a geomagnetic latitude of 
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+13.0° at a local time of 1200 on September 13, 1961, The geomagnetic field 
observed at this time on the Earth’s surface was moderately disturbed with 
several observatories reporting that a moderately severe magnetic storm 
commenced suddenly a few hours earlier at 1554 UT (Lincoln, 1962) . 

A noteworthy feature of the results shown in figure 42 is that the bound- 
ary of the magnetosphere is defined as much or more by the abrupt change in 
direction (note that ip changes by approximately 180°) as by the change in 
the intensity of the field. Just inside the magnetosphere, the magnetic field 
is relatively constant and of about twice the intensity of the extrapolated 
surface field, as assumed in the approximate calculation of the shape of the 
magnetosphere boundary outlined on figure 13. Beyond the boundary, the field 
is irregular, much as calculated by Auer, Hurwitz, and Kilb (1962) for 
Ma > 2. No evidence of passage through the bow wave is apparent in the data 
presented in figure 42, but it is quite likely that this is because the bow 
wave is farther from the Earth than the greatest distance for which data are 
indicated. 


The calculated results for Ma < 2 also have a counterpart in the data 
presented by Cahill and Amazeen (1963) and illustrated here in figure 43. On 
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Figure 43.- Explorer XII data resembling calculated results for antiparallel fields with M < 2 

(Cahill and Amazeen, 1963). A 

this occasion Explorer 12 crossed the magnetosphere boundary at a geomagnetic 
latitude of -18.2° at a local time of 0822. As in figure 42, there is again 
a sudden marked change in the direction of the magnetic field at the magneto- 
sphere boundary at 10.6 Earth radii. There is, however, no significant change 
in intensity. Farther from the Earth, the magnetic field appears to develop 
an oscillatory character that resembles, at least in broad outline, the calcu- 
lated results of Auer, Hurwitz, and Kilb (1961, 1962) for Ma < 2. In both 
figures 42 and 43, it may be observed that data from Explorer 12 more nearly 
resemble the calculated results for antiparallel rather than parallel fields. 
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Detailed analysis of the magnetometer data obtained over the entire use- 
ful life of Explorer 12, which extended from August 16 to December 6, 1961, 
has revealed many other important features of the outer magnetosphere and the 
surrounding flow. Mead and Cahill (1967) have examined the data covering the 
region from 3-1/2 to 13 Earth radii for seven passes during geomagnetical ly 
quiet times, and made comparisons both with the Jensen and Cain (1962) harmonic 
representation of the contributions of internal sources and the Mead (1964) 
bounded magnetosphere field model. Near the noon meridian and from 6 or 
7 Earth radii to the boundary of the magnetosphere, the magnitude and direc- 
tion of the observed field are shown to agree well with the results calcu- 
lated using the latter model. Near the dawn meridian, a significant twisting 
of the field lines out of the dipole magnetic meridian plane is observed, cor- 
responding to the tendency of the field to point in the antisolar direction in 
the southern hemisphere. On one pass north of the equator, the twisting is 
observed to be in the opposite sense from that expected, however. This result 
is interpreted as indicating that there may possibly exist large north-south 
magnetosphere asymmetries even when the solar wind is directed nearly 
perpendicular to the dipole axis. 

More than 70 distinct magnetosphere boundary crossings by Explorer 12 
have been analyzed extensively by Cahill and Patel (1967) . The location of 
the boundary was compared with the theoretical position determined by Mead 
(1964), displayed here in figure 19, normalized to a geocentric distance to 
the magnetosphere nose of 10.7 Earth radii. The measurements are generally 
consistent with the theoretical predictions, but the day-to-day boundary 
fluctuations and the lack of an effective plasma probe on Explorer 12 preclude 
making a precise evaluation of the theory. The thickness of the boundary was 
reported to be usually between about 20 and 300 km. Such distances are suf- 
ficiently small to justify the theoretical idealization of the boundary as a 
discontinuity surface for most purposes concerned with the external aerodynam- 
ics of the magnetosphere. On a few occasions, however, a sharply defined 
boundary is lacking, and thicknesses as great as several thousand kilometers 
are indicated if the boundary thickness is considered to be the entire dis- 
tance over which slow changes from magnetosphere to transition region values 
occur. 

Evidence of boundary motion was provided by the appearance of multiple 
boundary crossings, much like those observed earlier in the flight of Explorer 
10, in 22 of the 44 passages examined in detail. Motion of the boundary by 
more than 1000 km was inferred from measurements on 8 passes, but smaller 
motions (100 to 1000 km) appeared more frequently. Cahill and Patel (1967) 
have suggested that the boundary motion could, on occasion, be periodic. Evi- 
dence in support of the view that even the largest amplitude variations in the 
position of the magnetosphere boundary detected by Explorer 12 may be periodic 
is given by Freeman, Kavanagh, and Cahill (1967). The data are principally 
that obtained on the outward pass on August 12, 1961, when the satellite was 
in the fortuitous position to observe two complete cycles of a major boundary 
motion, together with more fragmentary results obtained on several other occa- 
sions. These data, limited though they may be, tend to suggest a periodicity 
of the order of 1/2 hour. Whether such a presumed periodicity represents the 
natural or resonant period of the magnetosphere, or the response to a driving 
periodicity in the solar wind remains uncertain. Quite independently of the 
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origin of these oscillations. Freeman, Kavanagh, and Cahill (1967) go on to 
suggest that they may be responsible for the 40 to 50 minute fluctuations found 
in the equatorial surface magnetic field by Pai and Sarabhai (1964) . 

On the other hand, absence of apparent boundary motion in half of the pas- 
sages examined indicates that gross boundary instability is not a characteris- 
tic of the forward part of the magnetosphere boundary. These results do not 
preclude such instabilities from occurring along the boundary of the 
magnetosphere tail, however. 

Cahill and Patel (1967) concluded from their analysis of the data that 
the boundary observations are generally consistent with the closed magneto- 
sphere model of Chapman and Ferraro described in this paper, although the 
boundary may, on occasion, be open in the manner suggested by Dungey (1963) . 
Additional support for these conclusions are provided in a related study of 
the magnitude of the magnetic field component normal to the magnetosphere’ 
boundary conducted by Sonnerup and Cahill (1967) . This field component was 
found, with a few notable exceptions, to be very small compared with the tan- 
gential component. In particular, they concluded there is little overall 
agreement with the open magnetosphere prediction of Levy, Petschek, and Siscoe 
(1964) that the magnitude of the normal component should be at least 10 per- 
cent and more likely 20 percent of the total field intensity. This result 
lends substantial observational support to the use of a hydromagnetic tangen- 
tial discontinuity surface, at which there is no normal component to the mag- 
netic field, to represent the steady-state magnetosphere boundary. Data 
consistent with a rotational discontinuity surface such as has been associated 
with the open magnetosphere model were observed on two occasions, both of 
which occurred while a geomagnetic storm was in progress. Under such circum- 
stances, the magnetosphere might be expected to be in a transient state of 
contraction or expansion, and the usual statements related to the steady-state 
magnetosphere boundary would not necessarily be expected to apply. 

Because apogee of Explorer 12 was only about 13 Earth radii, the bow wave 
was traversed on only three orbits at times when the magnetosphere was in a 
highly compressed state. Kaufmann (1967) has examined the data on each of 
these occasions. The location of the bow wave relative to the magnetosphere 
and the change in the magnitude of the magnetic field across the shock are 
concluded to be consistent with the results of fluid calculations such as 
those presented here in figures 22 through 29. The appearance of fluctuations 
in the magnetic field as the bow wave is crossed when traveling in the direc- 
tion toward the magnetosphere, and conversely, much as indicated by the small 
insert plots in figure 41, are interpreted as indicating the bow wave to be 
the source of disorder in the flow. The thickness of the bow wave was always 
thin compared with the overall dimensions of the magnetosphere, although dif- 
ficult to specify with precision in an unambiguous manner because of the com- 
plexity of its structure. The data indicate the time interval between the 
instant the magnetic field is first observed to deviate from its interplane- 
tary level and the instant when the quiescent field level first reaches the 
shocked value to vary between 30 seconds and the minimum resolving time of the 
measurements (less than 2 seconds) . The proper conversion of time intervals 
into distances is not known with certainty, but the latter may be estimated by 
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considering that the satellite is moving relative to the Earth at about 
1 km/sec, and the shock is thought to be moving past the satellite with an 
average speed of the order of 10 km/sec. Most of the data are judged to be 
consistent with the interpretation that this region consists of a shock much 
thinner than just indicated, plus several other waves. In addition, one or 
more large oscillations having a period of the order of several tens of sec- 
onds are sometimes observed near the shock wave. A number of possible explana- 
tions for these oscillations are discussed, but it is concluded that no 
decision can be reached as to which is correct. 


While the data from Explorer 12 together with the more fragmentary data 
from even earlier spacecraft served to suggest the possible existence of a bow 
wave upstream of the magnetosphere, it was the data from the magnetometers and 
plasma probes of Explorer 18, also known as IMP-1, that provided the definitive 
evidence. This spacecraft was launched from Cape Kennedy at 0230 UT on 
November 27, 1963, with an initial apogee of 31.02 Earth radii geocentric, 
perigee altitude of 192 km, and an orbital period of 3.934 days. The initial 
Sun- Earth- apogee angle was approximately 27° toward the dawn side of the Earth. 
The extremely eccentric orbit combined with the more than six months useful 
life of the spacecraft permitted mapping of the magnetosphere boundary and bow 
wave from near the subsolar point to far into the nighttime region of the 
Earth . 
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Clearly evident in much of the data are easily distinguishable differences 
between conditions in the magnetosphere, in the incident solar wind beyond the 
bow wave, and in the intervening transition region. As a result it became the 
practice of several of the experimenters to present summary plots showing the 
portion of the orbits in which the data displayed transition region character- 
istics together with the theoretical curves for the magnetosphere boundary and 
the bow wave reproduced here in figure 41. In the case of the magnetometer 
data, this region is recognizable by the presence of fluctuations more sub- 
stantial than in either the incident solar wind or in the magnetosphere, much 

as indicated by the results shown in 
figures 41 through 43. The differences 
are even more obvious in the data from 
the plasma probes . Figure 44 shows a 
sample of the latter, as presented by 
Wolfe, Silva, and Myers (1966), for the 
current output of their plasma probe 
during the inbound half of orbit 17. 

Only the data from the four lowest vol- 
tage or energy steps of the probe are 
plotted, since current was seldom 
observed in the energy windows above 
3740 volts (846 km/sec for protons) . 

This plot shows that no plasma flux was 
detected on this particular orbit when 
the spacecraft was nearer than about 
11.5 Earth radii. This is because the 
spacecraft was within the magnetosphere, 
and hence shielded from the flowing 
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Figure 44.- IMP-1 plasma probe output - orbit 17, 
January 31-February 1, 1964 (Wolfe, Silva, and 
Myers, 1966). 
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solar plasma. When it was beyond about 16 Earth radii, essentially all the 
plasma flux was observed by only one of the energy channels of the probe. 

This is because the spacecraft was in the solar wind, beyond the disturbing 
influence of the Earth, where the random thermal velocities of the particles 
were sufficiently small compared with the directed bulk velocity that the 
velocity distribution or energy spectrum of the particles was too narrow to 
bridge the gap between the energy windows of the plasma probe. The lack of 
current output in all energy channels when the spacecraft was near apogee is 
interpreted on the basis of related data for other times and of design features 
of the probe to be due not to a temporary absence of the solar wind, but to a 
shift in the bulk velocity to a value such that the sensible part of the par- 
ticle energy spectrum fell between the observing windows of the probe. 
Ground-based magnetometer data support this interpretation as well because 
they indicate that a gradual commencement geomagnetic storm occurred at approx- 
imately 0830 UT on January 31, 1964, when the spacecraft was at about 24 Earth 
radii. Presumably, then, the energy spectrum of the solar wind prior to this 
time would be such as to slip undetected between the 1700-volt and 600-volt 
channels of the instrument. 

Although the magnetosphere boundary and bow wave were not always so 
clearly defined in the data as in the case of orbit 17, it was usually possible 
to distinguish the data from the various regions. Figure 45 is a summary plot 

THEORY 



Figure 45.- IMP-1 plasma probe measurements of shock wave and magnetosphere boundary crossings. 

in which are drawn the portions of the first 29 orbits in which Wolfe, Silva, 
and Myers (1966) interpret the data from the plasma probe to indicate that the 
spacecraft is between the magnetosphere boundary and the bow wave. The dashed 
portions of some of the orbital segments represent possible positions of the 
boundaries in cases in which the termination of the transition region was 
uncertain. As discussed previously in the presentation of results from 
Explorer 12, such uncertainties could be due either to the gradual character 
of the change or the occurrence of multiple crossings resulting from the mag- 
netosphere boundary or shock wave moving back and forth past the spacecraft. 
This plot was derived by locating the transition region in terms of geocentric 
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distance and solar oriented ecliptic longitude. It therefore does not account 
for the latitude of the spacecraft during its traversal through this region. 

Superposed on the data are the theoretical positions for the magneto- 
sphere boundary and the bow wave calculated for flow with free-stream Mach 
number NL = 8 and ratio of specified heats y = 5/3. The bulk velocity Voo 
and density n^ of the solar wind are indicated to be 300 km/sec and 
10 protons/cm 3 , although any other combination of values having the same prod- 
uct n^v^ 2 would lead to identical curves for the magnetosphere boundary and 
the bow wave. Although the wide gaps between the energy channels and the 
limited directional information obtained by the plasma probe on IMP-1 prevent 
deduction of precise values for the density and bulk velocity of the solar 
wind, and hence the determination of detailed comparisons with theoretical 
results, the overall agreement of the gross features of the correspondence 
between the theoretical and observational results are clearly discernible in 
the results shown in figure 45. Although the criterion used to distinguish 
the transition region from the magnetosphere or the undisturbed solar wind are, 
of course, different, the same conclusion has been reached through examination 
of very similar appearing plots based on the magnetometer data from IMP-1 
(see, for example, Ness, Scearce, and Seek (1964) or Ness, Scearce, Seek, and 
Wilcox (1966) for extensive accounts) . 

Although perhaps not so central to the present discussion of the external 
aerodynamics of the magnetosphere, mention should be made of a number of 
related features of the solar wind and the magnetosphere revealed by the data 
from IMP-1. One of the most striking of these is the existence of a long mag- 
netic tail of the Earth and an enclosed magnetically neutral sheet extending 
far downstream from the Earth in the antisolar direction. Ness (1965) has 
reported that the diameter of the tail is approximately 40 Earth radii, and 
that magnetic field intensities of approximately 10 to 30 gammas are observed 
out to the satellite apogee of 31.7 Earth radii. The former is in deceptively 
good accordance with the theoretical results based on the Chapman-Ferraro 
theory, but the intensity and orientation of the magnetic field within the 
tail are very different from the current-free configuration visualized through- 
out the magnetosphere in that theory. The discovery of the neutral sheet and 
the attendant possibility of plasma instabilities have resulted in a vigorous 
development of theories to explain an extensive array of phenomena such as the 
aurora, gegenschein, day-night asymmetries of various sorts, and the formation 
of the Van Allen radiation belts. Accounts of many of these may be found in 
several other articles in this book. 

Since the magnetosphere tail displayed no appreciable attenuation out to 
apogee, data from IMP-1 could provide no insight into its ultimate length. 

More recently, however, Ness, Behannon, Cantarano, and Scearce (1967) have 
reported that both the tail and a relatively thin magnetic neutral sheet have 
been observed in the data from Explorer 33 satellite at geocentric distances 
as great as about 80 Earth radii. Behannon (1968) has also reported that data 
from the same satellite showed the Earth's bow wave to be still detectable at 
a geocentric distance of 75.7 Earth radii. The probable, although less cer- 
tain, detection of the geomagnetic tail at even greater distances ranging 
between 900 to 1050 Earth radii has been indicated by the plasma probe (Wolfe, 
Silva, McKibbin, and Mason, 1967) and magnetometer (Ness, Scearce, and 
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Cantarano, 1967) data from Pioneer 7. A coherent, well-ordered tail with 
embedded neutral sheet does not appear to have been observed, however, and it 
may be that the tail should be considered to be more like an irregular wake 
at this distance. Such an interpretation is consistent with the fact that 
measurements made by Mariner 4 spacecraft on a trajectory to Mars did not 
reveal any recognizable trace of the magnetosphere tail or wake as it passed 
through the anticipated region at a geocentric distance of about 3300 Earth 
radii (Van Allen, 1965) . 

Data from IMP-1 satellite were also valuable in confirming over an extend- 
ing interval of time the general Archimedean spiral structure of the average 
interplanetary magnetic field as indicated previously by solar wind theory and 
by more limited measurements made in space, particularly by Mariner 2 (Smith, 
1964). Particularly significant is the recognition by Wilcox and Ness (1965) 
of a long-lived magnetic sector structure corotating with the Sun in the inter- 
planetary medium. The boundaries of the sectors were frequently sharply 
defined in the magnetic records, but less so in the data from the plasma probes 
and other instruments. In general, enhanced values for the intensity of the 
interplanetary magnetic field, the solar-wind velocity, and geomagnetic activ- 
ity tend to be associated with the leading portion of each sector, and small 
values with the trailing portion. On the other hand, the density of the solar 
wind tends to be large in the leading and trailing portions of each sector, 
and small in the middle. The importance of these findings to the understanding 
of the long-recognized 27-day recurrence pattern of geomagnetic storms is 
immediately evident. 

More detailed information on the physical changes in the solar wind 
plasma as it crosses the Earth's bow shock wave has been presented recently by 
Argo, Asbridge, Bame, Hundhausen, and Strong (1967) using data from plasma 
probes on two Vela 3 satellites obtained on 13 shock crossings during July and 
August 1965. These satellites, which are in nearly circular orbits tilted at 
60° to the ecliptic plane, usually crossed the bow wave during this period at 
about 60° and 300° solar ecliptic longitude, and between ±56° ecliptic lati- 
tude, and at geocentric distances of approximately 18.5 Earth radii. In addi- 
tion to presenting the results for each crossing, the data were summarized by 
giving average values for the ratios of the velocities, temperatures, and den- 
sities on the two sides of the shock. They are v/v TO = 0.70, T/T^ = 24, 
p/Poo = 3.4, in which the subscript 00 refers to conditions in the incident 
solar wind and the symbol without subscript refers to conditions behind the 
bow wave. The corresponding theoretical values of Spreiter, Summers, and 
Alksne (1966) at the positions of the Vela shock crossings are quoted as 
v/Voo = 0.70, T/T^ = 12, p/p^ = 4, of which it is judged that only the tempera- 
ture fails to agree closely. Inspection of the latter reference or figure 23 
of this paper shows that this theoretical value for the temperature ratio is 
that associated with a free-stream Mach number M^ = 8 and ratio of specific 
heats y = 5/3, the combination usually selected for illustrative purposes, 
and that a more accurate value for the density ratio for this case is about 
3.8. The actual values for the speed and temperature of the incident solar 
wind listed for the 13 shock crossings correspond to substantially higher 
free-stream Mach number than 8, however. As a rough measure, it may be noted 
that the average of the Mach numbers for the 13 crossings is 11.3. If this 
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value is used for the Mach number together with the value 0.70 for v/v m , the 
expression for T/T^ shown on figure 23 leads, for y = 5/3, to T/T^ = 22.7. 
Alternatively, the individual observational values listed for the speed and 
temperature of the solar wind and the ratio v/v^ across the shock may be 
used to calculate a theoretical value for T/T^ for each crossing. Taking 
the average of these in the same way as the experimental values were averaged 
yields a temperature ratio of 23.8 instead of the originally quoted value of 
12. The new value is, of course, in virtually perfect agreement with the 
experimental value of 24, and tends to confirm the appropriateness of the 
value 5/3 for the ratio of specific heats as well as the usefulness of the 
aerodynamic model generally. 

The successful launching of Pioneer VI spacecraft equipped with a plasma 
probe having higher resolution than related instruments on earlier spacecraft, 
as well as a magnetometer and other instruments, provided an opportunity to 
make more detailed comparisons between the theoretical and observed character- 
istics of the solar wind and its interaction with the geomagnetic field than 
previously possible. Nature cooperated by providing a period of exceptionally 
steady flow of the solar wind during the time the spacecraft was traversing 
the region between the magnetosphere and the bow wave. On Earth, the geomag- 
netic field was also very steady with the planetary index Ap attaining a 
value of only 1 during the time of interest. The entire period surrounding 
the launch date was marked by low geomagnetic activity with the daily Kp 
sums on December 15, 16, and 17 being 2-, 1+, and 3-, respectively. 

Pioneer VI was launched on December 16, 1965, into an escape trajectory 
that remained very close to the plane of the ecliptic and appeared as illus- 
trated in figure 46 when viewed looking southward from over the north pole. 

The data from both the magnetometer (Ness, Scearce, and Cantarano, 1966) and 
the plasma probe (Wolfe, Silva, and McKibbin, 1968) show that the spacecraft 
crossed the magnetosphere boundary at 1256 UT at a geocentric distance of 
about 12.8 Earth radii and a solar-ecliptic longitude of 76.3°. It crossed 
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the bow wave somewhat more than four hours later at about 1712 UT when at a 
geocentric distance of approximately 20*5 Earth radii and solar-ecliptic 
longitude and latitude of 81.25° and -4.4°, respectively. 


Immediately after Pioneer VI passed the bow wave and emerged into the 
solar wind undisturbed by the presence of the Earth, the plasma probe detected 
a steady stream of particles coming from an apparent direction about 6° west 
of the direction to the Sun with a speed of about 280 km/sec and a number den- 
sity of about 11 protons/cm 3 (and presumably an equal number of electrons) . 
These values were used together with assumed values of 8 and 5/3 for the free- 
stream Mach number and the ratio of specific heats y to calculate the 

velocity and density of the flow and the intensity of the magnetic field all 
along the portion of the trajectory from inside the magnetosphere boundary to 
somewhat beyond the bow wave. The results are presented in figures 46, 47, 
and 48 together with those observed with Pioneer VI . The most apparent 
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Figure 47.- Comparison of Pioneer VI flow direction with calculated results for = 8, y = 5/3. 



Figure 48.- Comparison of Pioneer VI ion density and magnetic field intensity with calculated 

results for = 8, y = 5/3. 
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feature to be seen in the comparisons is the nearly perfect agreement between 
the calculated and measured locations of the magnetosphere boundary and the 
bow wave. It should be noted in this connection that although a particular 
value was selected for to make the calculations, the choice is not criti- 

cal. As described previously, the locations of the magnetosphere boundary and 
the bow wave, and also the contours of constant p/p^, v/v^, | B j / 1 | , and 
(T - T 00 )/v 00 2 , but not T/T^, tend to be virtually independent of for 

free-stream Mach numbers greater than about 5. 

The observational data for the peak flux ion velocity during the time 
Pioneer VI traversed the transition region, and for about a half-hour there- 
after, are presented in figure 46. In the incident solar wind where the ran- 
dom motions of the protons are small compared with their directed bulk 
velocity with respect to the Earth or the spacecraft, the peak flux ion veloc- 
ity corresponds very nearly to the bulk velocity of the flow, as can be seen 
by comparing the values for cm and w = v for a somewhat different case in 
figure 33. Behind the bow wave, however, substantial differences emerge. In 
order to illustrate their significance more clearly, two sets of theoretical 
curves are included on figure 46. The dashed line shows the bulk velocity cal- 
culated by direct application of the results for v/v^ presented on figure 23 
for a magnetosphere size determined in accordance with the results summarized 
in figure 15 for Voo = 280 km/sec and = 11 protons/cm 3 . The quantity 
represented by the data points is not the bulk velocity, however, but the 
velocity of the ions that produce the greatest current per unit energy incre- 
ment in the plasma probe. Of the various quantities calculated theoretically, 
this peak flux ion velocity is most nearly comparable to the speed c^ of the 
particles for which dn/n dC dQ is an absolute maximum. With the assumption 
that the solar wind is composed of fully ionized hydrogen, the relations given 
on figure 30 lead to the expression cjq = [w + (w 2 + 8kT/m) 1 ^ 2 ]/2 for c^ in 
terms of the bulkvelocity w and the temperature T. The solid line on fig- 
ure 46 indicates the values for the velocity cjq computed in this way for the 
conditions stated above. In confirmation of the more qualitative statement 
above, it may be seen that the values for w and cjq are nearly the same in 
the solar wind beyond the bow wave, and that substantial differences exist 
behind the shock wave where the temperatures are much higher and the bulk 
velocities somewhat lower than in the solar wind. In comparison, values for 
the peak flux ion velocity measured by Pioneer VI are slightly less than those 
calculated for cjq, although not so small as those calculated for the bulk 
velocity w. With the estimated 10-percent uncertainty in the determination 
of the experimental velocity distribution, however, the calculated results for 
both cm and w fall within the allowable range of values indicated by 
Pioneer VI. Aside from the matter of these modest differences in the absolute 
values, the general trend of slightly increasing velocity as the spacecraft 
moved along its trajectory from the magnetopause to the bow wave is 
represented very satisfactorily by the theoretical results. 

Figure 47 shows the corresponding results for the calculated and observed 
direction of flow. The agreement is generally quite satisfactory, although 
some significant deviations between the calculated and measured results may be 
seen in the data obtained near the magnetopause. In these plots, the flow 
directions are given in terms of a coordinate system fixed relative to the 
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Sun-Earth line as indicated in the small sketch. Because of the combined 
motion of the Earth and spacecraft, a flow that is radially outward from the 
Sun, as viewed in nonrotating heliocentric coordinates, would appear to 
approach the spacecraft with 0 = 0° and cp = 6° when = 280 km/sec. It may 

be seen that the flow actually observed by Pioneer VI when beyond the bow wave 

is nearly in accordance with these values. 

The corresponding comparison between the calculated and measured values 
for the density is shown in the top part of figure 48. Again the agreement is 
satisfactory . The plot in the lower part of the same figure shows the intensi- 
ties of the magnetic field measured by Pioneer VI (Ness, Scearce, and 
Cantarano, 1966) and those calculated using the theoretical results of figure 
24 together with the same values for M^, y , v M , and n M as described above, 
and a value of 3.5 gamma for selected to fit the measurements in the 

solar wind. Considering the complexity of the calculation of the magnetic 

field as one of the last quantities to be evaluated in a long chain of calcu- 
lations involving approximations at each step, the agreement is almost surpris- 
ingly good. It is, moreover, quite possible that some of the variations 
displayed by the observational data and not duplicated in the calculated 
results are the consequence of small changes in the intensity or direction of 
the interplanetary magnetic field, and could be satisfactorily reproduced by 
the theory if the interplanetary conditions were known. 

One feature of the results for which the latter statement is probably not 
true, however, is the pronounced peak or spike in the run of values for the 
velocity, density, magnetic intensity, and 0 component of the flow direction 
immediately downstream of the bow wave. The detection of such a localized 
irregularity is, moreover, not a unique property of the data from Pioneer VI, 
but a detail that has been repeatedly observed by other spacecraft in associa- 
tion both with the Earth f s bow wave and with propagating interplanetary shock 
waves. Although continuum gasdynamic or magnetohydrodynamic calculations by 
themselves provide no indication of such behavior, it appears likely that 
these spikes can be understood in terms of the calculated properties of strong 
collisionless shock waves such as illustrated in figure 41 for Alfven Mach 
numbers greater than 2. 

The preceding has not been intended to be an exhaustive or complete dis- 
cussion of the contributions of the various satellites, but rather an account 
of some of the findings that bear on the assumptions and predictions of the 
theories of the solar wind and its interaction with the Earth. In particular, 
data from several Vela satellites discussed elsewhere in this book and a num- 
ber of Soviet spacecraft have made important contributions, but space precludes 
further discussion here. In general, however, the results tend to be consis- 
tent with those described above, although certain features may be displayed 
with greater or lesser clarity depending upon the properties of the orbit, the 
type and resolution of the instruments, and conditions existing in space at 
the time of the flight. 
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SUMMARY OF STEADY ASPECTS OF THE INTERACTION OF THE 
SOLAR WIND AND THE GEOMAGNETIC FIELD 


Finally, in broad outline, one may summarize many features of the steady- 
state hydromagnetic interaction of the solar wind and the magnetosphere in 
terms of the sketch shown in figure 49. The solar wind approaches the Earth 
with both the Mach number and the Alfven Mach number much larger than unity. 



Figure 49.- The Earth as a magnetized obstacle in the solar wind. 

and with the ratio of the gas pressure p to the magnetic pressure B 2 /8tt 
of the order unity. The high Alfven Mach number assures that the equations of 
magnetohydrodynamics can be approximated with little error by the equations of 
gasdynamics for the flow and by the equations of frozen flux for the deforma- 
tion of the interplanetary magnetic field. The high Mach number indicates 
that the flow is similar to a supersonic or hypersonic flow and that a 
detached bow wave would be expected upstream of an obstacle. Behind the bow 
wave, the temperature and pressure are greatly enhanced over their values in 
the solar wind, while the magnetic field tends to be increased only moderately, 
and is even diminished in some regions. As a result, the gas pressure tends 
to be substantially greater than the magnetic pressure in the flow around the 
forward part of the magnetosphere. The Alfven Mach number tends to remain 
well above unity, although not so great as in the undisturbed solar wind. 

Inside the magnetosphere, the magnetic pressure is many orders of magni- 
tude greater than the gas pressure. The magnetosphere is bounded by a current 
sheath that adjusts its shape like a free surface so that the sum of the gas 
and magnetic pressure is the same on both sides at all points. Essentially, 
this is a balance between the gas pressure of the exterior flow and the mag- 
netic pressure of the confined and deformed geomagnetic field. With the 
introduction of the Newtonian approximation for the pressure of the flowing 
plasma, the determination of the shape of the magnetosphere boundary reduces 
to the classical Chapman- Ferraro problem based on the quite different concepts 
of particle, rather than fluid, motion. An incompatible condition is 
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disclosed in the vicinity of the neutral points, however, in which an ade- 
quacy of the Newtonian approximation allows the appearance of a sharply 
reflexed curvature in the shape of the magnetosphere boundary. Normally, a 
shock wave would be expected to form in the supersonic flow along such a sur- 
face. It is not possible for such a shock wave to extend to the magnetosphere 
boundary, however, because the associated jump in pressure cannot be matched 
by a jump in magnetic pressure in the stationary curl-free magnetic field 
characteristic of the magnetosphere. It is proposed, in view of the known 
behavior of fluid flows generally, that the flow separates from the boundary 
of the geomagnetic field and subsequently reattaches farther downstream. In 
this way there is formed an embedded region of hot, essentially stationary, 
plasma bounded partially by the free surface just described and partially by 
the cusped geometry of the magnetic field lines. Previous studies of the 
properties of similar cusped- geometry configurations considered for plasma 
containment for production of thermonuclear energy lead to the conclusion that 
a substantial flux of particles from this region leak into the magnetosphere 
from the vicinity of the end of the cusp. 

Although this paper is not concerned primarily with conditions in the 
interior of the magnetosphere, four prominent features of this region are 
included on this sketch for completeness. First of all, there is, of course, 
the approximately dipole geomagnetic field essential to the formation of the 
cavity, or in effect the obstacle, around which the solar wind must flow. 
Extending out to the magnetosphere boundary in the subsolar region, and to a 
lesser distance elsewhere, are the trapped charged particles of the extended 
Van Allen belt. Trailing far downstream from the Earth is the magnetosphere 
tail of great, although still undetermined, extent. Within this region the 
magnetic field is essentially rectilinear in form with the field vectors 
pointing away from the Earth in the southern half of the tail and toward the 
Earth in the northern half. The field diminishes in intensity and reverses in 
direction in a relatively thin and flat surface, the magnetic neutral sheet, 
that extends approximately down the center of the tail. The cross section of 
the tail thus resembles, at least for the first several tens of Earth radii, 
the Greek letter 6 with the magnetic field being approximately uniform in 
intensity throughout, but directed oppositely in the two halves. A somewhat 
greater than usual amount of plasma appears to be characteristic of the region 
near the magnetic neutral sheet, although the precise relation between the 
extensiveness of the plasma sheet and magnetic neutral sheet still remains to 
be determined. 


UNSTEADY ASPECTS 


Although all of the preceding discussion is concerned with steady-state 
phenomena, it has long been known that important transient effects exist. It 
was, in fact, observations of transient geomagnetic variations that led 
Chapman and others over 40 years ago into the pioneering studies that paved 
the way for the rapid understanding of many aspects of the data measured in 
recent years with spacecraft. Two distinctly different classes of major 
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irregularities in the interplanetary plasma flow are sketched on figure 50. 
One is a large plasma cloud of indefinite form and size expelled from the Sun 
by a flare or other explosive event. The existence of such a phenomena had 




GEOMAGNETIC 
FIELD CHANGES 

S.C. 

-avg-i — magnetic » 

storm 


time in days 



Figure 50.- Schematic drawings of two types of major irregularities in the solar wind and their 

geomagnetic consequences. 


been suspected long before the launching of the first satellite because of the 
frequent observation of a flare on the surface of the Sun followed a few days 
later by a set of geophysical phenomena described collectively by the term 
T, geomagnetic storm. " One characteristic feature of such an event is that dis- 
played in the upper right-hand part of figure 50 in which is depicted, in 
idealized form, the time variation of the horizontal component of the geomag- 
netic field observed in low and middle latitudes . The sudden commencement 
(s.c.) of many geomagnetic storms was originally thought to be the result of 
the impact of the plasma cloud on the geomagnetic field. More recently, how- 
ever, the tendency has been to associate it more with the impact of a possible 
bow wave traveling through the solar wind in advance of the plasma cloud. The 
second class of major irregularity in the solar wind is a long-lived rotating 
beam, possibly preceded by a bow wave. This type of irregularity, called an 
M-region beam, had also long been suspected because of a pronounced tendency 
for geomagnetic storms to recur in approximately 27-day intervals with each 
complete rotation of the more active regions of the Sun. Earlier discussions 
of both the transient plasma clouds and the long-lived rotating beams tended 
to regard them as intrusions into the otherwise effective vacuum of inter- 
planetary space. Since the announcement of the theoretical likelihood and the 
verification of the actual existence of the solar wind at all times, both the 
plasma clouds and the beams must be regarded as moving through a general back- 
ground of presumably less dense and slower moving flow of solar plasma. 

Although spacecraft have now provided definitive evidence of the actual 
existence of such irregularities in the solar wind, their probable existence 
could have been inferred long ago from photographs of comets that occasionally 
display substantial changes in the direction of the tail in the course of a 
few hours. A famous pair of photographs illustrating such a change in the 
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appearance of Comet Morehouse during a three-hour period on September 30, 1908, 
is shown in figure 51. Clearly visible is the substantial change in the direc- 
tion of the tail much as though a front across which there is a sudden shift 
in the wind direction were moving along the length of the tail. 



Figure 51.- Yerkes Observatory photographs of Comet 1908 III, Morehouse, showing change in 

3 hours on September 30, 1908. 

Impact of an irregularity of any sort in the solar wind upon the magneto- 
sphere may be expected to initiate a variety of transient phenomena. Analysis 
of the steady- state theory indicates that a change in the momentum flux 
PooVqo 2 of the solar wind may be expected to result ultimately in a change in 
the size of the magnetosphere, and probably in the establishment of oscilla- 
tions that would gradually damp with time. It is reasonably well ascertained 
observationally , although the precise details of the connecting mechanisms 
remain obscure, that such changes in the solar wind are also the cause, 
directly or indirectly, of several other important geophysical phenomena. 

Major fluctuations of the contents of the Van Allen belts may occur, accom- 
panied by significant heating and enhanced ionization of the upper atmosphere 
sufficient, in more intense events, to impair radio communications, particu- 
larly in the polar regions. Auroral phenomena also appear to be affected sig- 
nificantly with increased activity and movement to lower latitudes being 
associated with certain phases of a geomagnetic storm. 

In contrast to Chapman and Ferraro's concentration on transient phenom- 
ena, most of the more recent analyses of the interaction of solar plasma and 
the geomagnetic field, particularly those in which fluid models are invoked, 
are concerned with the steady-state problem. As a result, most of the quanti- 
tative theoretical analyses of large-scale transient features of the interac- 
tion are based on the older particle models that do not include the presence 
of a bow wave upstream of the magnetosphere. To the extent that disturbances 
do indeed propagate through the bow wave without excessive qualitative change, 
and the Newtonian formula, when expressed in a reference frame fixed with 
respect to an element of the magnetosphere boundary, continues to provide a 
reasonable approximation for the pressure, these results may still provide a 
guide to many features of the interaction. 
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The rapidity with which the magnetosphere can respond to a major change 
in the solar wind is a matter of interest in many applications. Although the 
qualitative manner in which the geomagnetic field would carve a cavity in a 
finite cloud of solar plasma advancing through a vacuum toward the Earth was 
one of the first topics discussed by Chapman and Ferraro (1931) , the actual 
quantitative solution of their problem was not given until recently. In fig- 
ures 52 and 53 are shown the results calculated numerically by Spreiter and 
Summers (1965) for two cases in which an extensive cloud of plasma having an 
initial density of 5 protons/cm 3 (and an equal number of electrons) and veloc- 
ity of 500 km/sec advances toward and engulfs the Earth. In both cases, the 
plasma cloud is assumed to have a flat face when at a geocentric distance of 
20 Earth radii. In figure 52, this face is normal to the direction of motion; 
in figure 53, the plasma face is inclined 30° thereto. The solid lines show 
the shape of the face of the plasma cloud in the equatorial plane at various 
times. The dashed lines show the final steady-state configuration of the 
boundary defined by the approximate solution described and illustrated in 
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Figure 52.- Successive positions in equatorial plane of front of advancing plasma cloud. 
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Figure 53.- Successive positions in equatorial plane 
of oblique front of advancing plasma cloud. 


figures 13 through 16. The results 
are free of surprises and quite con- 
sistent with the early qualitative 
description of Chapman and Ferraro. 
Quantitative features of particular 
interest are the rapidity with which 
the boundary approaches its steady- 
state form as the. plasma front sweeps 
by, the smallness of effects of obliq 
uity of the plasma front, and the 
tendency to form waves that drift 
downstream along the boundary as it 
nears the steady-state configuration. 

The results presented in dimen- 
sional form in figures 52 and 53 are 
convenient for providing a quick idea 


58 



of the quantitative features of the interaction for a representative set of 
parameters, but lack generality as presented. The corresponding results for 
plasma clouds of speed and density different from those indicated may be 
determined quickly, however, by application of simple scaling laws inherent in 
the basic equations. They indicate that one must only alter the distance and 
time scales in accordance with the relations rb/ r a = (Pst^/Pst^) 1/6 and 
tb/t a = (Pst a /Pstb^ 1/6 C v . a / V *b) in which Pst = PooVoo 2 is the stagnation 
pressure, subscript a refers to the values used to calculate figures 52 and 
53, and subscript b refers to any other set of values. 

The preceding problem is important because of its relative simplicity, 
conceptual value, and historical interest. The continual existence of the 
solar wind limits its actual usefulness, however. We consider next that a 
steady solar wind has existed for a sufficiently long time (of the order of 
10 minutes or more on the basis of the preceding results) that the boundary 
separating the geomagnetic field and the solar plasma is well established in 
the steady-state configuration associated with a value Psti ^ or stagna- 
tion pressure. We consider further that a discontinuity in momentum flux 
advances through or with the solar wind plasma and that the discontinuity sur- 
face is planar over lateral distances larger than the magnetosphere before it 
begins to interact with the magnetic field of the Earth. Behind the discon- 
tinuity surface, conditions are considered to be steady and uniform, but are 
characterized by a larger stagnation pressure Pst 2 • Ultimately, therefore, 
the magnetosphere boundary will assume a configuration geometrically similar 
to that possessed originally, but smaller in the ratio (p s ti/Pst 2 ) 1/6 - The 
time history for the transition between the two steady-state configurations 
has been calculated by Spreiter and Summers (1965) and the results are shown 
in figure 54 for a case in which the initial steady state is the same as the 
final steady state arrived at in the case portrayed in figures 52 and 53. The 
particular discontinuity considered is one that drifts along with the local 
velocity of the solar wind and across which there is a fourfold change in the 
density. The initial and final steady-state solutions are indicated by the 



time, sec 

Figure 54.*- Successive positions in equatorial plane of oblique density discontinuity in 

solar wind. 
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outer and inner dashed lines. As in the previous examples, the solid lines 
represent the shape of the discontinuity surface at the stated times. The mag 
netosphere boundary at each of these times is given in part by the portion of 
the solid line that is interior to the outer dashed line. The remainder of 
the magnetosphere boundary is still stationary in this approximation in the 
initial position indicated by the outer dashed curve. As in the previous 
examples, the results show that the magnetosphere boundary collapses rapidly 
to its final configuration. Although several minutes are required for the dis 
continuity surface to sweep past the forward part of the magnetosphere, each 
element of the boundary essentially completes its movement from initial to 
final position in less than a minute. 

It is not difficult to visualize the fluid counterparts of the problems 
described above, or to formulate their mathematical representation using the 
equations of magnetohydrodynamics. Because of the formidable difficulties 
encountered in the solution of these equations, however, quantitative discus- 
sion of such problems has hardly begun. Dryer, Merritt, and Aronson (1967) 
have described certain aspects of the dynamic interaction between the magneto- 
sphere and bow wave system and a blast wave propagating through the solar wind 
with the aid of a laboratory experiment in which a blast wave produced by a 
shock tube is directed downstream through an operating supersonic wind tunnel 
containing a simplified model of the magnetosphere in the test section. 

More localized aspects of the interaction such as the initial changes in 
a gasdynamic or magnetohydrodynamic shock wave or other discontinuity upon 
impact with the bow wave can be treated immediately by standard methods 
described in texts such as Courant and Friedrichs (1948) , Landau and Lifshitz 
(1959), Oswatitsch (1956), and Jeffrey and Taniuti (1964), While methods, and 
even solutions in some cases, are immediately available from these and other 
sources, their application to the magnetosphere remains largely unexploited. 

As an example of the ease of such an application, we may consider the amplifi- 
cation of a plane sound wave upon head-on impact with a strong shock wave. 

Such an idealization might be considered sufficient to make the first esti- 
mates of what happens to a small disturbance in the solar wind upon encounter- 
ing the forward part of the Earth's bow wave. For this problem, Landau and 
Lifshitz (1959) show that the shock wave is set into oscillation, and that 
both a sound wave and a contact surface wave are to be found in the region 
downstream from the shock wave. Also shown is that the ratio of the amplitude 
5p2 of the pressure oscillations in the transmitted sound wave to the ampli- 
tude 6pi of the incident sound is proportional to the pressure ratio across 
the stationary shock wave. More precisely, 

6p^/6pi = {y + [2y(y - l)] 1/2 }" 1 (p 2 /Pl) 

or approximately 0.317p 2 /pi for y = 5/3. For a strong shock wave, 

P 2 /P 1 = 2yMi 2 /(y + 1) , or 5M} 2 /4 for y = 5/3, in which M} should be equated 
to the incident Mach number based on the component of velocity normal to the 
shock wave. With Mx = 8, as might occur at the nose of the Earth's bow wave, 
these expressions indicate that p 2 /px = 80 and 6p^/6px = 25.4. This result 
illustrates the immense amplification that small disturbances may experience 
as they propagate across the Earth's bow wave. It also furnishes a possible 
alternative explanation to that of collisionless shock wave theory illustrated 
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in the small inserts of figure 41 for the increased irregularities in the data 
from plasma probes and magnetometers behind the bow wave compared with corre- 
sponding data obtained in the incident solar wind. The different speeds of 
propagation of the transmitted sound wave and the contact surface wave pro- 
duced by interaction of a sound wave and the bow wave offers the possibility 
of further richness of consequences because of interactions that may occur 
between the two transmitted wave systems. Burgers (1946) has shown, for 
example, that there will be a simple wave train of velocity and pressure dis- 
turbances propagating at the local speed of sound, and also a simple wave 
train of entropy disturbances carried along with the local velocity of the 
flow. The interaction of the density and temperature oscillations associated 
with both of these wave trains leads, in addition, to the development of oscil- 
lations with new frequencies by a process similar to that by which beats are 
produced. Details of the interaction of the bow wave and the wide variety of 
waves and other irregularities possible in the solar wind plasma remain 
largely undetermined, as do the characteristics of their propagation through 
the flow field exterior to the magnetosphere, and beyond to the interior of 
the magnetosphere and ultimately to the surface of the Earth. 

The structure of the bow wave remains a topic surrounded by considerable 
uncertainties. Attention has been given in this text to the possibility that 
the structure may be understood in terms of the collisionless shock wave cal- 
culations of Auer, Hurwitz, and Kilb (1961, 1962) and Rossow (1965, 1967). 
Another possibility has been developed extensively recently by Tidman (1967a, 
b) in which a turbulent shock wave is considered to be formed through the 
action of ion waves produced as a result of two-stream instabilities of inter- 
penetrating ion streams. In spite of these and several other proposals, it 
does not seem unfair to conclude that the precise determination of the 
detailed mechanism by which the bow wave is formed in the collisionless plasma 
of the solar wind remains largely unknown at the present time. 

The stability of the magnetosphere boundary is another topic that has not 
been touched upon in the preceding discussion. It is apparent that the bound- 
ary must be stable in at least some gross sense, since otherwise the generally 
good agreement between observations and the steady-state theory could hardly 
be expected to occur. Numerous theoretical studies of the question of the 
stability of the magnetosphere boundary have been published over the last 10 
years or more, but the results are often contradictory. Dungey (1958, 1963) 
considered that the flow of solar plasma along the magnetosphere boundary 
would generate surface waves in the same way that the wind generates waves on 
water through the action of a Kelvin-Helmholt z instability. He and Parker 
(1958) gave a theory of such waves, and concluded that the magnetosphere 
boundary is unstable. The applicability of the results to the magnetosphere 
is not assured immediately, however, because of the neglect of known features 
of the phenomena such as the compressibility of the plasma, the curvature of 
the boundary, and the effects of nonlinear terms in the governing equations. 
Spreiter and Summers (1965) proceeded differently by basing their analysis on 
the strict application of the Chapman- Ferraro theory. They found that, if the 
wave length and amplitude are sufficiently small that curvature and second- 
order effects can be disregarded, all perturbations, except those having wave 
fronts alined with the direction of the local magnetic field within the mag- 
netosphere, damp exponentially with time and drift along the magnetosphere 
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boundary with the tangential component of the solar wind velocity. Alined 
waves, which neither damp nor amplify in this approximation, are examined fur- 
ther by inclusion of curvature and higher order effects. Analysis shows that 
curvature introduces a destabilizing effect in small regions of opposite to 
normal magnetosphere boundary curvature in the vicinity of the neutral points 
and a stabilizing effect elsewhere. With the elimination of the regions of 
reverse curvature near the neutral points, as described in connection with 
figure 39, even this exceptional region of instability disappears. Buffalano 
(1967) has considered this same model further, and shown that convection of a 
boundary perturbation along the magnetosphere boundary acts also to limit the 
growth of any instability. On the other hand, several others including Talwar 
(1964), Sen (1965), Lerche (1966) have developed analyses based on fluid con- 
cepts, and shown, for assumed conditions that may or may not adequately repre- 
sent the actual conditions, that the magnetosphere boundary is unstable. More 
extensive accounts of this important and controversial subject have been given 
recently by Atkinson and Watanabe (1966) and Buffalano (1967) . 


CONCLUDING REMARKS 


Three broad conclusions emerge from this review of the external aerody- 
namics of the magnetosphere. One is that substantial progress is being made 
in the development of methods for the calculation of the properties of both 
the solar wind and its interaction with the geomagnetic field. Another is 
that calculations based on fluid concepts of magnetohydrodynamics and gasdynam- 
ics are consistently of almost surprising accuracy considering both the ques- 
tions surrounding the use of continuum methods and the introduction of 
simplifying assumptions necessary to achieve solutions at this time. Finally, 
one cannot avoid the conclusion that a considerable amount of further analysis 
is required to bring a corresponding clarification to many problems that still 
remain unsolved. 

In the solar wind, the true nature of the heating mechanism in the corona 
that drives the solar plasma away from the Sun remains unknown in detail. A 
discussion of the manner by which the solar wind terminates at great distances 
from the Sun and merges with the interstellar medium, perhaps through the 
action of one or more terminal shock waves, has only begun. The ultimate con- 
sequences of such inquiry may be of considerable importance in understanding 
the interaction of stars and their gaseous environment in space. The quanti- 
tative analysis of the properties and propagation of transient effects is also 
of great interest and must be understood in much greater detail than at pres- 
ent for an adequate comprehension of a vast array of phenomena that come under 
the general heading of solar- terrestrial relations. 

Any list of unsolved problems in the interaction of the solar wind and 
the geomagnetic field must include the quantitative description of the mag- 
netosphere tail, and the development of a complete solution of the shape of 
the magnetosphere boundary and the properties of the surrounding flow in a 
completely self-consistent manner. It does not appear likely, however, that 
the solution of the latter very difficult problem would yield results substan- 
tially different from those that can now be calculated. An important problem 
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in need of clarification is the question of the stability of the magnetosphere 
boundary. The structure of the bow wave and the mechanism by which it creates 
or amplifies irregularities in the flow is a topic of great intrinsic interest 
and important work needs to be done to provide even a rudimentary understand- 
ing. The properties and propagation of transient effects all the way from the 
solar wind, or even more ambitiously, from the surface of the Sun, to the 
Earth remains a subject of vital importance that is still relatively undevel- 
oped. Also of interest, although barely discussed as yet, is the extension 
and adaptation of the understanding achieved for the Earth to other objects in 
the solar system. This may possibly be straightforward in the case of a mag- 
netized planet such as Jupiter, or very different in the case of essentially 
unmagnetized objects such as the Moon, Mars, and Venus. 

Finally, perhaps the most intriguing question of all relates to the very 
success of the continuum models in dealing with the. most rarefied gas pres- 
ently accessible to direct measurement. It is "Why does collisionless, 
slightly magnetized plasma act so much like an ordinary collision-dominated 
gas?" It is not that one is entirely devoid of explanations, but those 
usually offered are somewhat unsatisfying and superficial involving little 
more than the simple replacement of the mean free path by some shorter inter- 
action distance such as the gyro radius. The proper resolution of this ques- 
tion must certainly be ranked among the more important goals to be sought if 
the environment provided by interplanetary space is to live up to its 
frequently touted potential as a laboratory for the study of rarefied plasmas. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, October 27, 1967 
129-02-08-03-00-21 
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